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The  purpose  of  this  report  is  to  examine  the  problem 
of  scattering  of  electromagnetic  waves  by  an  open,  perfectly 
conducting  surface.  Specifically,  to  formulate  the  problem 
as  a boundary  value  problem,  to  convert  the  latter  into  a 
problem  in  integral  equations,  and  to  show  that  it  can  have  at 
most  one  solution.  In  trying  to  reach  these  objectives  we  were 
guided  by  the  methods  employed  to  achieve  the  same  ends  for  the 
problem  of  scattering  by  a closed  surface.  In  this  latter 
case  there  are  two  items  of  special  concern:  the  types  of 

surfaces  and  the  types  of  linear  current  densities  induced  on 
these  surfaces  that  would  allow  for  a successful  completion 
of  the  aforementioned  tasks.  In  the  case  of  an  open  surface 
these  same  items  appear  but  in  a more  complicated  form:  in 

the  interior  of  the  surface  things  are  not  any  different  from 
a closed  surface;  near  the  edge,  however,  additional  require- 
ments must  be  imposed  both  on  the  surface  and  on  the  induced 
linear  current  density.  Though  it  is  not  immediately  apparent 
from  reading  this  report,  much  of  the  time  devoted  to  its 
preparation  was  spent  in  finding  the  right  class  of  open 
surfaces  and  current  densities  for  which  the  objectives  could 
be  accomplished. 

For  convenience  the  report  has  been  divided  into  two 
parts.  In  Part  II  we  deal  with  problem  described  above.  In 


surface,  of  the  quasi-stationary  form  of  the  integral 
representations  for  the  scattered  fields  obtained  in  Part  II. 
This  part  is  essential  in  completing  the  proof  that  the 
boundary  value  problem  is  equivalent  to  a problem  in  integral 
equations . The  apparent  reversal  in  the  order  of  the  two 
parts  is  intentional,  for  it  is  in  what  we  call  Part  I that 
all  the  restrictions  on  the  surface  and  some  on  the  current 
density  come  into  play.  The  other  possible  arrangement, 
i.e.  making  Part  I an  appendix  of  the  main  problem,  would 
necessitate  bringing  in  these  restrictions  at  the  beginning  of 
the  report  without  any  possible  explanation.  Again  for 
convenience,  we  have  supplied  each  part  with  its  own  title, 
abstract,  and  introduction. 
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PART  I I 

li 

H 

The  Behavior  of  Potentials  with  Singular  Densities  [| 

Near  the  Edge  of  an  Open  Surface. 


Abstract 

t! 

This  work  deals  with  the  behavior  of  a simple-layer 
potential  and  its  spatial  derivatives  near  the  edge  of  an 
open  surface.  Types  of  open  surfaces  and  density  distributions, 
singular  near  the  boundary  of  such  surfaces,  are  determined 
for  which  these  potentials  satisfy  a finiteness  of  energy 
condition.  The  resulting  estimates  are  useful  in  establishing 
integral  representations  for  the  electromagnetic  fields 
scattered  by  a perfectly  conducting  open  surface. 


A. 


Introduction. 


It  has  long  been  recognized  (Jones,  1964;  Soitunerfeld, 

1964)  that  the  conditions  required  for  scattering  problems 
involving  closed  surfaces  to  be  well-posed  are  not  sufficient 
in  case  the  surfaces  are  open.  In  fact,  it  has  been  shown 
(Jones,  1964)  that  in  order  that  the  open  surface  scattering 
problem  have  at  most  one  solution  (if  at  all) , then  one  more 
condition  is  required  of  the  scattered  fields,  be  they  acoustic 
or  electromagnetic.  This  condition  requires  that  the  fields 
behave  in  a specified  way  in  the  vicinity  of  the  edge  of  the 
open  surface  and  is  usually  referred  to  as  an  edge  condition. 
Though  a condition  of  this  type  is  dictated  purely  from  the 
mathematics  of  the  problem,  it  has  also  a physical  meaning: 

It  is  intimately  connected  to  what  is  known  as  the  finiteness 
of  energy  condition  which  requires  that,  in  a volume  devoid 
of  sources,  the  energy  content  should  go  to  zero  with  the 
volume. 

In  this  work  we  will  examine  whether  a simple-layer 
potential  and  its  first  derivatives  satisfy  such  an  edge  or 
energy  condition.  This  question  came  about  from  studying  the 
problem  of  finding  integral  representations  for  the  electro- 
magnetic fields  scattered  by  a perfectly  conducting  open 
surface.  Once  found,  these  representations,  which  are  in 
effect  the  Stratton-Chu  formulas,  must  be  tested  as  to  whether 
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they  satisfy  an  edge  condition.  It  can  be  shown  that  this  can 
be  accomplished  by  examining  a simple-layer  potential  whose 
density,  defined  on  the  open  surface,  is  allowed  to  increase 
beyond  bound  (though  in  a specified  way)  in  the  approach  to 
the  boundary  of  the  surface.  Specifically,  we  will  show  that 

I the  line  integral  of  the  square  of  the  potential,  and  also  of 

its  spatial  derivatives,  around  a circle  centered  on  the 
I boundary  of  the  open  surface  and  lying  on  a plane  perpendicular 

to  the  boundary  exhibit  a certain  order  behavior  with  respect 
• to  the  radius  of  the  circle.  In  Section  B we  will  introduce  the 

j types  of  open  surfaces  to  be  considered.  In  Section  C we  will 

describe  the  class  of  allowable  densities  and  will  state  the 
problem  and  the  results.  In  Sections  D through  G we  will  prove 
the  results  stated  in  Section  C,  while  in  Section  H we  will 
offer  some  concluding  remarks.  Some  detailed  computations  will 
be  left  for  Appendices  A through  D. 

Before  closing  we  mention  that  problems  of  this  type  have 
' i been  studied  in  great  detail  in  two  dimensions  for  Cauchy-type 

t 

4 

' integrals  whose  density  is  defined  on  an  open  curve 

(Muskhelishvili , 1953;  Gakhov,  1966).  As  above,  the  density  is 
allowed  to  grow  beyond  bound  near  the  end-points  of  the  curve. 

I Using  these  results,  Hayashi  (1973)  was  able  to  obtain  order 

relations  near  the  end-points  of  the  curve  for  the  scattered 
fields  for  the  two-dimensional  Dirichlet  problem  for  the 
Helmholtz  equation.  Some  criticism  of  his  results  is  offered 


1 


in  Section  H. 


a unit  normal  on  z = F(x,y) , we  choose  the  one  whose 
z-component  points  in  the  same  direction  as  the  positive 
z-axis . 

3 

Definition  2 . A surface  S in  R is  a smooth  open  surface 
if 

(a)  S is  a smooth  surface 

(b)  for  any  two  points  in  which  do  not  belong  to 

S there  exists  at  least  one  continuous  curve 
connecting  them  which  does  not  have  any  points  in 
common  with  S. 

(c)  S has  as  its  boundary  a simple  closed  and 
rectifiable  curve  C. 

Definition  3.  A smooth  surface  S is  an  orientable  surface 
if  it  is  possible  to  cover  it  with  a family  of  coordinate 
neighborhoods  as  in  (1)  in  such  a way  that  if  a point  M of 
S belongs  to  two  neighborhoods  of  this  family,  then  the  change 
of  coordinates  has  positive  Jacobian  at  M. 

This  definition  says  in  effect  that,  at  each  point  in  the 
common  part  of  two  overlapping  neighborhoods,  the  normals 
for  each  neighborhood  point  in  the  same  direction. 

Definition  4.  A surface  S in  R^  is  a bounded  surface  if 

there  exists  a sphere  that  contains  S in  its  interior. 

3 

Definition  5.  A surface  S in  R 


is  a connected  surface 
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if  any  two  of  its  points  can  be  connected  by  a continuous 
curve  on  S. 

Definition  6.  A surface  S in  R is  a simply  connected 
surface  if  it  is  connected  and  if  every  simple  closed  curve 
on  S can  be  continuously  deformed  into  a point  of  S. 

The  type  of  surface  we  are  interested  in  is  a bounded, 
simply  connected,  orientable,  smooth  open  surface  with  some 
additional  conditions  on  its  boundary  as  well  as  the  portion 
of  the  surface  near  the  boundary.  The  condition  of  boundedness 
together  with  that  of  smoothness  result  in  the  surface  having 
area  (cf.Carmo,  1976).  The  condition  of  connectedness  is  not 
really  necessary.  The  actual  surface  may  be  composed  of  a 
(finite)  number  of  disjoint  surfaces.  Each  such  surface  again 
could  be  multiply  connected  provided  that  the  "holes"  in  it  are 

3 

not  composed  of  a single  point  in  R , and  that  the  curves 
bounding  these  "holes"  as  well  as  the  surface  near  them 
obey  the  conditions  described  below.  As  it  will  become  clear 
later  on,  however,  an  extension  of  the  ensuing  results  to  such 
surfaces  is  immediate,  and  for  this  reason,  as  well  as  to  avoid 
cumbersome  notations,  we  restrict  ourselves  to  a simply  connected 
surface . 

We  now  turn  to  the  description  of  the  boundary  as  well  as 
of  the  surface  near  it.  From  this  point  on  we  call 
the  surface  S,  its  boundary  C,  and  its  closure  S (S  = SUC) . 
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The  properties  of  the  previous  paragraph  apply  now  to  S 
and  not  just  S.  Thus,  if  M is  a point  of  C,  we  have 
the  representation  (1)  for  the  portion  of  S inside  the 
sphere  of  radius  d with  center  at  M.  We  require  next  that 

at  each  point  of  C there  exists  a well-defined  tangent 

line  which  we  identify  with  the  x-axis  of  the  local  coordinate 
system.  We  take  the  x-axis  to  be  positively  oriented  with 
respect  to  the  z-axis,  so  that  the  y-axis  is  pointed  towards 
S.  With  respect  to  this  coordinate  system,  we  impose  the 
following  condition  on  S: 

For  all  points  M of  C there  exist  fixed  positive 
numbers  A,  a',  and  b'  such  that 

|F(x,y)-F(0,y)  | < a|x|^'*'^,  |x|  < a'  < d,  0 < y < b'  < d (2) 

where  y ^ y > 0 and  may  depend  on  the  point  M of  C. 

This  condition  can  be  given  a geometric  interpretation. 

The  unit  normal  vector  for  the  portion  of  the  surface  under 
consideration  is  given  by 

^ -F  x-F  y+z 
X Y 

/l+F  ^+F  ^ ■ . 

X y 

Since  from  (2)  F^(0,y)  =0,  0 £ Y £ < d,  we  have  that  on 

the  yz-plane,  the  unit  normal  is 


I 


I 

1 


\ 

I 


[ 


i 
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-F^(0,y)y+g 

n = - - 

/l+Fy^(0,y) 

But  this  is  precisely  the  unit  normal  vector  to  the  two- 
dimensional  curve  z = F(0,y),  0_<y£b'  <d,  provided  its 

arclength  is  increasing  with  y.  Thus,  Eq.  (2)  requires 
that  the  three-dimensional  curve  z = F(0,y)  behaves  as  a 
two-dimensional  curve  in  0£y<b'<d,  which  in  turn  says 
that  the  curve  has  torsion  zero,  (or,  equivalently,  its 
osculating  plane  is  the  yz-plane) . This  lack  of  torsion  can 
be  loosely  restated  in  terms  of  the  surface  itself  as  follows: 

At  X = n,  the  surface  is  not  allowed  to  twist  about  the 
curve  z = F(0,y),  0_<y_<b'<d. 

With  respect  to  the  boundary  C of  S,  we  require  that 
it  is  twice  continuously  differentiable  with  respect  to  its 
arclength  and  that  the  portion  of  C intercepted  by  the  sphere 
of  radius  d and  center  a point  M of  C has  a projection 
on  the  xy-plane  described  by  the  pair  of  equations. 

X = f(s),  y = g(s),  -Sq  £ s £ s^^  (3) 

where  s represents  arclength  measured  from  the  origin  and 
increasing  in  the  direction  of  the  positive  x-axis.  The  functions 
f and  g possess  continuous  second  derivatives,  bounded  third 
ones,  and  satisfy  the  conditions. 
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f(0)  = g(0)  = 0,  f'(0)  = 1,  g'(0)  = 0,  f ' (s)  ^+g ' (s)  ^ = 1.  (4) 

We  note  here  that  the  curve  described  in  (3)  does  not  intersect 
itself  since  the  point  (x,y)  is  in  the  domain  of  F. 

We  call  the  conditions  described  by  Eqs . (2),  (3),  and 

(4)  the  open  surface  edge  conditions  and  we  make  the  following 
definition. 

Definition  7.  A surface  S in  is  a regular  open  surface 

if  it  is  a bounded,  simply  connected,  orientable,  smooth  open 
surface,  and  if  it  satisfies  the  open  surface  edge  conditions 
(2)  - (4) . 

This  is  the  type  of  surface  we  will  be  considering  from  now  on. 

C . Description  of  the  problem. 

The  objective  here  is  to  study  the  behavior  of  a simple- 
layer  potential  and  its  first  derivatives  near  the  boundary  of 
a regular  open  surface.  The  potential  is  of  the  form 


U(P) 


; 

■'s 


h(M) 

R(P,M) 


dS 


(5) 


where  R(P,M)  is  the  Euclidean  distance  between  the  points  P 
and  M,  M e S,  P ^ S,  and  h is  the  density  function  defined 
on  S.  As  is  well  known  this  type  of  potential  has  been  studied 
in  great  detail  for  S a closed  Lyapunov  surface,  and  h 
Holder-continuous  on  it  (Gunter,  1967;  Kellogg,  1953)  . In 
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the  present  case  the  surface  is  open  and  the  density  will  be 
allowed  to  grow  beyond  bound  in  the  approach  to  the  boundary. 

In  order  to  define  the  density  function  precisely,  we  need  to 
introduce  a new  curvilinear  system  of  coordinates. 

With  center  a point  of  C we  erect  rectangular 

coordinates,  as  described  in  the  previous  section,  and  we  focus 
on  the  portion  of  A'  contained  in  the  rectangle  |x(  £ a', 

0 £ y £ b' • The  unit  tangent  vector  at  a point  of  the  plane 
curve  (3)  is 

= f ' (s)x  + g'  (s)y,  (6) 


while  the  unit  normal  vector  is 


riQ  = zx£q  = -g'(s)x  + f'(s)y. 

The  position  vector  then  to  a point  (x,y)  of  A'  is 
xf(s)  + yg(s)  + nQ(s)p,  p £ 0 


(7)  . 


(8) 


from  which  we  get  the  transformation 


X = f(s)-g'(s)p,  y = g(s)+f'(s)p 
which  will  be  one-to-one  provided  the  Jacobian 


9 (x,y)  ^ 
9 (s , p) 


K (s) 


(9) 


(10) 


is  not  zero.  Here,  <(s)  is  the  curvature  of  the  curve  in 
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(3) , i.e. 

K(s)  = f ' (s)g"(s)  - g' (s)f"(s) . (11) 

By  the  assumption  on  f and  g,  k is  bounded  and  we  can 
choose  p sufficiently  small  so  that  the  Jacobian  is 
strictly  positive.  We  assume  that  this  is  the  case  for  the 
curvilinear  rectangle 

A={(s,p)  : |s|_<a,  0£p£b}c[-a',a']  x [0,b']>  (12) 

The  density  function  h is  now  defined  as  follows: 
it  is  a real-valued  function  defined  over  S and  satisfying 
the  following  two  conditions 

a)  In  every  closed  and  connected  subset  S'  of  S, 
bounded  away  from  C,  the  function  h is  continuous. 

b)  In  A of  (12),  i.e.,  near  and  on  the  boundary,  the 
function  h is  of  the  form 

h(M)  = , 0<ct<l,  a + y>l,  (x,y)  e A (13) 

P 

where  u is  the  index  introduced  in  (2)  and  a is  a 
function  defined  on  A and  satisfying  the  Lipschitz 
condition 

|o(M^)  - o(M2)|  < a|Mj^-M2|,  M^,  M2  g a . (14) 

Definition  8.  A real-valued  ^unction  h defined  over  a 
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regular  open  surface  S is  said  to  be  a regular  density 
function  for  the  surface  if  it  satisfies  conditions  (a)  and 
(b)  above. 

As  mentioned  at  the  beginning  of  this  section,  the 
point  P in  (5)  is  near  the  boundary  C of  S and  is  chosen 
as  follows:  Given  a point  Mq  of  C,  we  make  it  the  origin 

of  the  rectangular  coordinate  system  described  above.  We  take 
the  point  P to  be  a point  on  the  yz-plane  with  coordinates 
(y',z')  and  such  that  0 < (y'2+z'2)^^^  = p'  < b.  With  this 
point  we  also  associate  polar  coordinates  with 


y'  = p'  cos  4>' , z'  = p'  sin  4> ' , p'  > 0,  0 £ <l> ' £ 2tt.  (15) 

In  the  succeeding  sections  we  will  examine  the  behavior  of 

the  line  integral  of  the  square  of  U(P)  as  well  as  the  square 

2 2 2 

of  its  spatial  derivatives  on  the  circle  y'  +z'  = p'  in  the 

limit  as  p ' ■+  0 . The  results  can  be  summarized  in  the  following 
two  theorems . 

Theorem  1 : Let  S = S U C be  a regular  open  surface  (Def.7  ) 

and  h a regular  density  function  (Def.  8)  defined  on  S. 

Let  U(P)  be  the  simple-layer  potential  defined  in  (5),  where 
P = (0,y',z'),  M = (x,y,z)  with  the  coordinate  system  as 
in  the  preceeding  paragraph.  Let  also 


U^(P) 


xh  (M) 

3 


S R(P,M) 


dS. 
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Then, 


U(P)  = 0(1),  Uj^(P)  = 0(1),  p'  0,  0 < <p'  < 2n. 


Theorem  2:  Let  S and  h be  as  in  Theorem  1 and  define 


UjtP)  . 


ds,  u,(p) . 

S R(P,M)  ^ 


( 2-z ' ) h (M) 
S R(P,M)^ 


dS, 


Then , 


■2ir 


p'  lUj  (P)  rd<{) 


^d(h'  = i 


0(p'^“‘^)  , 


0 < a < 1 


P'  ^ 0 , 


m2, 


[0 ( p ' (log  p ' ) ) , a = 1 


where  j = 2 , 3 . 


We  remarlc  here  that  the  estimates  in  Theorem  1 are  stronger 
than  those  of  Theorem  2.  Several  attempts  were  made  to  obtain 
estimates  for  U2  and  li)ce  those  for  U and  U^,  but 

were  not  successful.  In  each  case  the  bounds  would  depend  on 
the  angle  (fi ' on  a way  that  would  make  them  non-square- 
integrable  with  respect  to  ij) ' . 

D.  The  behavior  of  the  simple-layer 


We  split  the  integral  in  (5)  into  two  integrals  by  writing 


U(P) 


S(A) 


h(M) 

R(P,M) 


dS  + 


S-S (A) 


h(M) 

R(P,M) 


dS, 


(16) 
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where  S(A)  is  that  neighborhood  of  S about  whose 

projection  on  the  xy-plane  is  the  region  js]  £a,  0 < p_<b. 

The  second  of  these  integrals  is  continuous  in  a neighborhood 
of  the  point  Mq  , and  its  limit  as  p ' 0 is  equal  to  the 

analagous  integral  obtained  when  P is  replaced  by  Mq.  The 
first  integral  can  be  written  as 


R(P,M) 


dS  = 


/l+F  ^+F  ^ dx  dy 

A p^-'^  (P,M)  ^ y 


— A+F  (l-fc(s)p)dsdp, 

A p “r(P,M)  ^ ^ 


where  R(P,M)  = (y-y ' ) ( z-z ' ) ^ and  a*(s,p)  = 

o(x(s,p),y(s,p)).  We  then  have 


R(P,M) 


dS  < M' 


■ , , r'  = /x^+(y-y')^  (18) 

P r ' 


where  M'  = max{ |a*(s,p) (l-K(s)p) ( /l+F^^+F^^  : (s,p)  e A}.  If 

in  (18)  we  expand  x and  y about  s = 0,  we  find  that 
(see  Appendix  A) , 


4^+(p-y')^ 


d+0(e')  ) , e’  ^ 0 , 


provided  that  b £ £'/6(|k|  +1)  <1,  a < e'/2(l+Mj^)  < 1, 


<|b  < 1,  with  K = K(0)  , = 6M+M'^,  and  where  M is 
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1 


the  maximum  of  |f"|'  I?"!*  can  then 

write 


S(A) 


h(M) 

R(P,M) 


dS 


< M(1+0{e')) 


p dsdp 


^ /s^+(p-  |y'  1 ) ^ 


e'  0 


(20) 


This  last  integral  exists  and  is  of  0(1)  as  y*  -*■  0,  (see 
Appendix  C) . We  thus  have, 

U(P)  = 0(1),  p'  0,  0 <(})'£  27t  (21) 

and,  hence, 

f2TT 

p'  |U(P)  I ^d(|)'  = 0 (P-  ) , p'  0.  (22) 

^0 


E.  The  behavior  of  Uj^(P). 

In  this  section  we  deal  with  the  integral 


Uj^(P) 


xh  (M) 
R(P,M) ^ 


(23) 


As  with  eq.  (16),  this  integral  can  be  split  into  two  integrals, 
the  second  of  which  is  well  behaved.  For  the  first  one  we  write. 


xh  (M) 

S(A)  R(P,M)^ 


dS 


[a(x,y)sec  ip  - o(0,|y'|)sec  ' ] 

J A 


a-1 

— ^ dxdy  + 
R(P,M) 


+ a (0 , I y ' I ) sec  ip  ' 


xp 


a-1 


A R(P,M) 


j dxdy. 


(24) 


I 
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where, 

sec  [p  = »^+F  ^+F  ^ |(x,y),  sec  ip'  = /l+F  ^+F  ^ |(0>ly'|)'  (25] 

X y ^ y 

i.e.  \p  is  the  angle  the  z-axis  makes  with  the  normal  to  the 
surface  at  (x,y)  while  ip'  is  the  corresponding  angle  with 
the  normal  at  (0,|y'l).  Since  F is  twice  continuously 
differentiable  we  have  that 

/"2  2 

I sec  Ip  - sec  t]; ' | £ const . /x“^+ (y- | y ' | ) ; (26) 

similarly,  since  o satisfies  (14) 


|o(x,y)  - a(0,(y'()|  £ const,  + (y-  y’ I ■ 


Combining  the  two  equations , we  obtain 


1 a (x,y)  sec  ip  - a (0  , | y ' | ) sec  ij; ' | _<  c /x^+  (y-  | y ' | ) ^ , (28 ) 


c a constant.  We  then  have  that  the  first  integral  on  the 
right-hand  side  of  (24)  is,  in  absolute  value,  less  or  equal 


f l.xlp“'Vx^(y-|^’ll^  ^ 

U (x2+(y-y-)^)^/^  JA  /2.,~ 


^ /x^+(y-|y'  1)^ 


This  last  integral  is  of  the  type  encountered  in  (20)  and  is 
of  0(1)  as  y ' -*■  0 . 


r 
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For  the  remaining  integral  in  (24)  we  write 


xp 


a-1 


A R(P,M) 
where , 


-j  dxdy  = 


xp 


a-1 


R(P,M)  R 


dxdy  + 


xp 


a-1 


A R' 


dxdy,  (30) 


R 


= /x^+  (v-y ' ) 


(y-y')  +(F(0,y)-z')^ 
For  the  expression  in  the  brackets  we  write 


1 

1 1 

' 

1 1 1 

R(P,M)  ^ 

R^ 

,R(P,M)  ^ 

2 ^2 
R{P,M)  R(P,M)R  R J 

For  the  first  term  of  this  expression  we  have 


R^-R(P,M) ^ 


and 


so  that 


R(P,M)  R R(P,M)R[R(P,M)+R]  ' 

R^  - R(P,M)^  = [F(0,y)-z']^  - [F(x,y)-z’]^ 
= [F(0,y)-F(x,y) ] [F(0,y)-z'+F(x,y)-z'j 


1r^  - R(P,M)^|  < |F(x,y)-F(0,y) 1 (R+R(P,M) ) 


and , by  ( 33) , 


R(P,M) 


|F(x,y) -F(0,y) | 
R(P,M) R 


(31) 


(32) 


(33) 


(34) 


1 

r 


I 


Combining  this  result  with  (2),  we  have  for  (32) 
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1+U 


R(P,M)  R 


3|  - R(P,M)R 


R(P,M) 


2 R(P,M)R 


R 


3A  X 


.4 


( 


where  r'  is  given  by  (19).  The  first  integral  on  the  right- 
hand  side  of  (30)  is,  in  absolute  value,  less  than  or  equal  to 


3A 


I 2+y  a-1 

Lx  I P 


,.4 


dxdy  < 3A 


p dxdy 


+ (y-|y'  I ) 


2-m 


(36) 


The  last  integral  exists  and  is  bounded  for  all  values  of  y' 
provided  a + y > 1,  (see  Appendix  C) . 


In  order  to  evaluate  the  remaining  integral  in  (30)  we 
resort  to  the  coordinates  (s,p)  introduced  in  (9).  We  thus 
write , 


xp 


A R 


j-  dxdy  = 


f p (f  .(.s)  -g'  (s)  p)  (i_K(s)  p)dsdp.  (37) 


R' 


Expanding  the  integrand  about  s = 0 (see  Appendix  B) , we  have 


xp 


a-1 


A R 


j-  dxdy  = (1+0 (e ' ) ) 


p“  ^ [ (1-Kp) s + i 5s^] 


— j 5 ^ (l-<(s)  p)dsdp,  e'-»-0  . 

A [s  +(p-y)  ')^+(F(0,p)-z')^]-^/"^ 


(38) 


From  (11) 


ic(s)  = <(0)  + [f ' (s.)g'"(s.)-g' (s.)f  "'(s.)]s,  0 < s.  < s. 


so  that 


As  with  (16)  , this  integral  can  be  split  into  two  integrals, 
the  second  of  which  is  well  behaved.  For  the  first  one  we 
write , 


^S(A) 


(y-y* )h(M) 
R(P,M)  ^ 


[o(x,y)sec  \p  - a(0,|y'|)sec  ' ] 

J A 


/ I X Ol-l 

dxdy  + 

R(P,M) 


+ 0 (0  , I y ' I ) sec  ' 

+ CT  (0  , |y ' I ) sec  ip  ’ 


(y-y' ) p“  ^ 

1 1 

A 

R(P,M)  R _ 

dxdy  + 


, , . a-1 

dxdy. 

A R"^ 


(44) 


The  first  two  integrals  on  the  righthand  side  can  be  treated 
in  the  same  way  as  in  the  previous  section.  For  the  last 
integral  we  have,  according  to  Appendix  B and  the  last  section. 


IizyLLp 

A R^ 


a-1 


= (l+0(e')) 


f'" 

0 • 

dxdy  = 

a-1 


(P-V’ ) (l-<p) p 


-a  [s^+(p-y')^+(F(0,p)-z')^] 


-2^2  + 0(1)  = 


=(l+0(e') )(2a 


(p-v') (l-Kp)p 


a-1 


° [ (p-y'  ) ^+(F(0,p)-z')  /a^+(p-y'  ) ^+(F(0,p)  -z’  ) ^ 


dp+0 (1) 


If  we  let  M*=max{ |F(0,p)|=0£p<b},  we  then  have 

(p-y')^  + (F(0,p)-z')^  < (b+|y'|)^  + (M*+|z'|)^- 


(45) 


From  this  inequality  together  with  the  requirement  that 


b+)y'  I 


< ^ , 
/2  /2 


(46) 


we  have  that 


(p-y')^  +(F(0,p)-z')^  < a^ 


(47) 


which  allows  us  to  expand  the  radical  in  (45)  into  a power  series 
to  obtain 


(y-y’ ) P 


a-1 


dxdy 


(l+0(e') ) 


(p-y') (1-Kp)p^~^ 
(p-y')^+(F(0,p)-z')^ 


dp+0 (1) 


(48) 


We  observe  that  the  last  integral  in  (48)  can  be  written  as 


b , I ^ , n ~ V 

■ <P-Y'I  ^ do  - Re 

0 (P-Y')  +(F(0,p)-z')^ 


b a-1,,  ~ . 
p (l-<p) 

T-w 


dp 


(49) 


where , 


T = p + iF(0,p),  0 < p < b,  and  w = y'+iz' 


(50) 


To  study  the  complex  integral  in  (49)  we  introduce  polar 
coordinates 

p + iF(0,p)  = re^®,  0 = 6(r)  (51) 


with 


= /p^+F(0,p)^  > p, 


(52) 


From  Appendix  D we  also  have 
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dr  > ^ dp , r < — p . 
2 " /3 


Taking  into  consideration  that  l^jp  < If  we  have  that 


b 01  - 1 , , , 

P (1-kp) 


b a-1 


'b  „a-l 


< :(1H-K|P).  dp  < _8  ' - r-  -dr 

■I  0 Ire^^-p'e^*^  | /3  ■'0  [ re^^-p ' 


where  rj^  = r(b) 


From  the  discussion  up  to  now  we  see  that 


2 [f^ti  r“"^dr 

|V-(P)P  < C.  ■■  . , 

2 ' — 1 „ I 10  ,i(I)'| 

■'0  re  -p'e  ^ 


+ 0(1) 


0 re  -p ' e 


0(1),  p'  ^ 0, 


where  is  a constant.  Since  the  aim  is  to  estimate  the 

integral 


p'  Iv  (P)  I ^d(J)' 


we  start  by  estimating  the  expression 


I(p’)  = 


2tt  r f b 


a-1  , 
r dr 


^ /r^  + p'^-2rp'  cos  ( (J) ' -0  (r ) ) 


d0’f  (57) 


which  we  can  write  as 


■f 
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I(P')  = 


2t7  r^b 

P'd4.' 

0 j 0 J 0 


(rq)°‘  ^drdc 


^r^+p  ' ^-2rp  ' cos  ( 4> ' -0  (r)  ) /q^  + Q ' ^-2qp  ' cos  ( 4) ' -0  (q)  ) 


We  consider  the  integration  on  the  (r,q,(J)')  space.  The  integrand 
is  not  defined  on  the  surface  obtained  by  revolving  the  curve 
r = p',  0 = (|) ' about  the  (|)'-axis;  otherwise  it  is  continuous 

on  (0,rj^]  X (0,rj^]  x [0,2ti].  Hence,  the  integrand  is  measurable 
(cf . Sobolev,  1964) , We  show  below  that  the  iterated  integral 
with  the  order  of  integration  reversed  exists.  By  the  Tonelli- 
Hobson  theorem  (cf.  Apostol,  1974),  it  is  then  equal  to  the 
original  integral.  We  have  that 


I(p')  < P' 


^b  f’^b 


0 J 0 


(rq)“  ^drdq 


d^T r _dJJ 

0 r^+p'^-2rp'  cos  (4i  ’ -0  (r)  ) 0 q^+p'^-2rp'  cos  ( <{■. ' -0  (q)  ) 


where,  above,  we  used  the  Cauchy-Schwarz  inequality.  If  we  make 
the  substitution  y = ~ 6 we  see  that  the  integrands  in  the 

last  two  integrals  are  periodic  in  y with  period  2n . We  can 


then  write, 


f b c b 1 

I ( p ' ) _<  P ' (rq)  drdq 

Jo  Jo 


0 r^+p'^-2rp'  cos  <)) ' J 0 q^+p'^-2qp'  cos  ' 


-23- 


The  integrals  with  respect  to  4>'  can  be  evaluated 
by  using  the  following  integral  representation  for  Legendre 
functions  of  the  second  kind  (Magnus  et  al.,  1966,  p.  186) 


Q**,  (z)  = i(2ir)  ^(z^-1)**  [ (z-cos  t)  ^dt, 

J n 


Since  (Magnus  et  al.,  1966,  p.  172) 


2,.  -hi 

(z)  = i(j)^(z  -1) 


we  have  that 


(z-cos  t)  ^dt  = ■tt(z^-I)  ^ 


From  (63)  we  can  compute  the  integrals  in  (60) : 


2 2 

0 r +p'  -2rp'  cos  4)' 


0 r^+p'^ 

-JFF^  - cos  4)’ 


Equation  (60)  then  becomes 


I ( p ' ) < 2irp  ' 


(rq)°^"^drdc 


0^0  /|2  ,2i|2  ,2 

/|r  -p'  ^ \q  -p’ 


^ r^'^dr 
0 lr^-p'^1 
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We  make  the  change  of  variables  r = p'£  to  obtain  for  (65) 


2a-l  f ^dr 

I(p’)  < 2ttp'^°‘  ^ 


In  evaluating  (66)  we  consider  the  cases  0 < a < 1 
and  ct  = 1 separately.  For  0 < a < 1, 


b/p-  ^ p c«-^dc 

0 k^-il^  " •'o 


and  the  last  integral  exists  as  an  improper  integral.  In  fact, 
through  simple  transformations  it  can  be  shown  that 


•oo  a— 1 

g d^ 

0 


a/2-1, 


4 t Ml-t)"^"  "dt  + 4 t ^(1-t)  dt 


"a/2-^3 


_ iLfl  a]  „[l  l-a]‘ 

■ 2L®[2'  2J^®[2'  ^JJ  ' 


where  B stands  for  the  Beta  function.  We  then  have  that 


x/»i\  i2a-ir„fl  al,„fl  l-a)l  „ „ ^ , 

I(p')  ± 2 P ® 2'  7 2'  ~2~  ' 0 < a < 1, 


For  the  case  a = 1 we  compute 


^b/p' 


2 I 5^ 

r-ii  ^ 


i"  ii 


■b/p' 


IT  , T b 

^ + log  ■— V- 

2 ’ p ' 


r-1 


< J + 1 + log  ^ ; 


T 
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so  -that  by  substitution  into  (66)  we  obtain 


I ( p ' ) < 2iTp  ' 


1 . IT  , T b 

1 + _ 4-  log  p- 


a = 1 


:71) 


To  estimate  the  second  term  in  (56) , i.e.  the  one 
resulting  from  the  second  term  in  (55)  , we  use  the  Caucliy- 
Schwarz  inequality  and  the  results  for  I(p');  thus, 


■2tt 

d(|)'p' 

0 


a-1  j 

r dr 


2rp ' cos 


< /2ttp'I(p')  < 


ri 

a 

M 

1 

P 

B 

l2' 

2 

+B 

12'  2 J 

< ^ 


2irp  ’ 


1 + ^ + log 


lb 
P’ J 


p ' , 0 < a < 1 

a = 1. 


(72) 


Finally , 
f 2it 


p’  |V2(P)  I d(^'  = 


0 (p 


I 2ci  1 % 


■ 

f \ 

r. 

2‘ 

0 

P' 

log 

D 

Ip'i 

0 < a < 1 


, a = 1 


0.  (73) 


G.  The  behavior  of  U2(P). 

In  this  section  we  examine  the  integral 

U (P)  = - dS.  (74) 

^ JS  R(P,M) 

As  with  (16),  this  integral  can  be  split  into  two  integrals, 
the  second  of  which  is  well-behaved.  For  the  first  one  we 


write , 


^ '>  ds=  [a  (x,y)  sec  -a  (0  , |y ' I ) sec  <p']  ^ , dxdy  + 

S(A)  R(P,M)  Ja  R(P,M)-^ 


+CJ  (0  , ly ' I ) sec  i(^ ' I p°‘  ^ (z-z  ' ) i - ~ dxdy+  [ ^ dxdy 

|Ja  R(P,M)  R Ja  R 


The  first  two  integrals  on  the  right-hand  side  can  be  handled 
as  in  the  previous  two  sections.  For  the  third  one  we  write 


dxdy  = f dxdy  + f dxdy. 

A R^  Ja  p-^  “r  JA  p-^  “r-" 


By  (2),  the  first  of  the  integrals  on  the  right  hand  side  is, 
in  absolute  value,  less  than  or  equal  to 


a-1  I I 1+U 


■X  P X P dxc 

i L_j dxdy  < A 

J A R ^ A f / 2^,  I , 

[/x  +(y-|y' 


d“l  j j 

p dxd^ 


(77) 


The  last  integral  is  of  the  same  type  as  the  one  in  (36).  For 
the  last  integral  in  (76)  we  write,  using  Appendix  B, 


F(0  y)-z'  f p“"hF(0,p)-z'+  i 

■rl^  3-  <3rfy  = U+0(E))  — ^ 5 ^ 2-372 

A p-^  “ R^  Ja  [s^+(p-y')^+(F(0,p)-z')^]'’/^ 


(l-K(s) p)dsdp  = (1+0  (G) ) 


[F(0,p)-2‘  ] (l-<p)  p^ 


-a  [s^+(p-y') ^+(F(0,p)-z') 


dsdp+0  (1)  ). 
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As  with  (45) -(48),  we  have 


F(0 ,y)  -2  * 

. l-a„3 

■'  A p R 


dxdy 


(l+0(e) )]2| 

1 ^ 

'0  (p-y')  +(F(0,p)-z')^ 

We  observe  now  that 


(79) 


fb 
• 0 


[F(0,p)-z']  (1-Kp)p“  ^ 

(p-y')^+(F(0,p)-z')^ 


dp  = -Im 


fb 
•'  0 


p°^  ^(l-<p) 
T-W 


dp , 


(80) 


and,  hence,  all  the  results  of  the  previous  section  can  be 
used  here.  The  sought-after  bound  is  the  same  as  in  (73),  i.e. 


2tt  2 

p'  IV  (P)  1 d<))'  = ^ 
0 ’ 


0(p'^“  ^) , 0 < a < 1 


2^  p' 

, a = 1 


O p log 


P ' 


0 . 


(81) 


Conclusion . 


The  main  results  of  this  work  are  Theorems  1 and  2 of 
Section  C.  In  formulating  these  theorems  at  the  outset  we 
were  certain  as  to  the  class  of  density  functions  we  wanted 
to  have  included  (Def.  8)  but  not  as  to  the  class  of  open 
surfaces.  Condition  (2)  on  the  surface  as  well  as  the  condition 
a + M > 1 in  (14)  were  not  anticipated  and  were  made  necessary 
by  the  manner  we  proved  these  theorems.  It  will  be  interesting 
to  try  to  construct  alternate  proofs  which  do  not  employ  these 
two  conditions.  Another  point  for  investigation,  which  we  also 
brought  up  in  Section  C,  is  whether  order  relations  with 
respect  to  p'  can  be  ol.tained  for  U2  and  which  would 


( 

I 


, I 
' 1 
r 1 

I 


f 


1 


) 


i 
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lead  to  the  same  result  as  in  Theorem  2.  With  respect  to 
this  last  point,  we  wish  to  mention  that  the  proof  Hayashi 
(1973)  gave  for  the  tv;o-dimensional  case  is  not  satisfactory. 

The  reasons  for  it  are  outlined  below. 

Consider  a simple,  open,  and  bounded  arc  L on  the 
yz-plane  and  make  one  of  the  end  points  the  center  of  a coor- 
dinate system  as  in  (15)  with  (say)  the  y'-axis  tangent  to 
L at  the  end-point.  Consider  the  function  u(p ',({)')  which  is 
the  combination  of  a simple-  and  a double-layer  potential  with 
densities  defined  on  L.  Hayashi  shows  that  u = 0(1) 
as  p ' 0 . For  3u/9p'  he,  however,  reasons  as  follows: 

"If  we  assume  that  3u/3p'  = 0(p’'^)  , p ' ->■  0 , then  a > -1 
is  necessary  in  order  that  u is  bounded  when  p'  -*•  0." 

This  statement  is  not  wrong  but  on  the  other  hand  it  is  not  a 
proof  that  3u/8p'  = 0(p'°''),  p'  -*-0,  as  claimed  in  his 
Theorem  A. 2.  What  has  been  proven  is  that  if  u(p')  =0(1) 
and  3u/9p'  = 0(p'^^)  , then  a > -1.  We  wish  to  make  the  point 
here  that  even  though  his  conclusion  is  incorrect,  this  in  no  way 
affects  the  results  of  his  paper  for,  as  in  our  case,  he  is 
interested  in  an  integral  of  the  type 

■ 2tt 

p ' I u ( p ' , ij) ' ) I d(J) ' 

0 

vanishing  in  the  limit  as  p ' 0 and,  we  believe,  though 


i 


1 


I 

I 

1: 
■ 1 


I 

i 


I 


we  have  not  proven,  that  this  is  indeed  the  case. 
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L.  Proof  of  Equation  (20)  . 


We  start  with  Equation  (9)  which  we  repeat  here  for 


convenience : 


X = f(s)-g'(s)  , y = g(s)+f'(s)  . 


(A.l) 


Without  loss  of  generality  let  s > 0.  Expanding  x and  y 
about  s = 0,  we  find 


X = (1-k(0)p)s  + i[f  " (Sj^) -g"  ' (Sj^)  p]  s^  , 


y = P + i-[g"  (s2)+f"  (S2)p]s^, 


(A. 2) 


where  0 < s^  < s,  0 < S2  < s.  We  then  have 


2 22  2~~222 
: +(y-y')  = s^+(p-y')^+(-2K  +K^^p^)s^  + 


3 (p-y' ) s^+ (1-Kp)  as^  + |-(a^+3^)s^, 


where  k = k(0) , and 


(A. 3) 


a = a(s,,p)  = f " (s, ) -g" ' (s, ) p,  6 = 6(s_,p)  = g" (s~ ) +f " ' (s_ , p) p , 


(A. 4) 


2 2 

Letting  B stand  for  s +(p-y')  while  A for  the  rest  of 
the  right-hand  side  of  (A. 3),  we  can  write 


^ £ 1 -2  Kp+ K^p^  1 + I 3 1 I s 1 + 1 1- Kp  I I a I I s^  + |-(a^+S'^)  s'^ . (A.  5) 


l,r.2.;:2,_2 


I 


r 


'I 


1 1 
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Assuming  that  f and  g have  second  and  third  derivatives 
bounded  by  the  constant  M,  we  get 


a|  £ (1+p)  , I 6 1 £ M(l+p) , 


(A. 6) 


so  that 


£ |-2+kp| |k|p+M(1+p) |s|+|1-kp|M(1+p) [si  + y M^(l+p)^|sl^ 


_<  (2+kp)|k|p+2M(1  + p)|s|  + |ic|pM(1+p)|s|  +iM^(l  + p)^|s|^ 


(A. 7) 


Further  requiring  that  p<l,  |k|p<1,  |s|<1,  we  get 

for  (A. 7) 


< 3 I k| p + 4M| s| +2M| s| +M^ I s I = 3 | k | p + 6M | s | +M^ | s | = 3 | k | p+M^ | s | , 

(A.  8) 


where  = 6M+M  . Letting  0 < e < 1 be  given,  we  have  that 

[A/b]  < e provided  that 


I 


6(|k|+1) 


< 1,  I s 1 < 


2 (1+M^) 


< 1,  k|p  < 1.  (A. 9) 


We  now  examine  the  expression 


/B+A  /b  A+A/B  /B 


1 + 


/1+A/B 


- 1 


(A. 10) 


For  0 < A/B  < E we  have 


while  for  -e  < — < 0 , 
B — 


1 

/l+A/B 


1 


< 


1_ 

/in 


1 . 


since,  for 


< e 


lim  — 
£-0+  ^ 


1 

/l+A/B 


(A. 11) 


we  have  from  (A. 10) 

— — = — [l+0(e)],  £ 0+  . (A. 12) 

/B+A  /B 

Without  loss  of  generality  we  assume  that  the  conditions 
in  (A. 9)  for  (A. 12)  to  hold  are  satisfied  in  the  curvilinear 
rectangle  A of  (12) . 


1.  The  expression  for  1/R'^  in  terms  of  the 


coordinates  (s,p). 


We  wish  to  rewrite  1/R  of  (31)  in  terms  of  the  coordinates 
(s , p) . Let 


u = F(0,y)  = F(0,g(s)+pf ' (s) ) 


(B.l) 


Since 


= g'  (O)-^pf " (0)  = 0 


we  have 


2 

u = F(0,p)  + j = F(0,p)  + j ys^,  0 < < s 

^ 


where 


^ p,  = 

3s  s=s. 


= i_Fl0xX)_  (g-  (s3)+pf"  (S3) 

3y  s=S3 


(B.2) 


3F(0,y) 

3y 


(g- (s3)+pf" • (S3)  ) 


(B.3) 


From  (B.2)  we  have 


[F(0,y)-z'l^  = [F(0,p)-z'l^  + y[F(0,p)-z']s^  + j , 


so  that 


+ (y-y')^  + (F(0,y)-z')^  = B'  + A',  (B.4) 


where 


B'  = B + (F(0,p)-z']^ 


A'  = A + y[F(0,p)-z']s^  + i 


(B.5) 


(B.6) 


and  A and  B have  the  meaning  given  to  them  in  Appendix  A. 
Then, 


<_  1 -2kp+kp^  I + ( 1 3 I + 


s I + 1 1-K  p 1 1 a 1 I s I + j(a^+B^+Y^) 


and  by  (A. 5)  and  (A. 8) 


jr  <3|k|p+M^|s|+ 


I ^ 1 ~2  2 ,r.  ns 

s|  + J y s . (B.7) 


From  (B.3)  with  all  derivatives  bounded  by  M, 


1y1  1 (1  + P)  (1  + P)  < (M'^+M'^)  (1  + p) < (M+l)^(l+p)^  < 4(M+1) 


.3,„2 


3,,  ,..2 


(B.  8) 


so  that  for  (B.7)  we  have 


gT-  < 3|k|p  + [M^  + 4 (M+1)  ^ + 4 (M+1)  1 s 1 = 3|<lp  + | s j . 

(B.9) 

Choosing  p and  s according  to  (A. 9)  with  substituted 

by  M2 / we  have  that  |A'/B'|  < e. 
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Appendix  C.  Computation  of  the  integrals  in  Equations  (20) 
and  ( 36 ) . 


The  integral  in  (20)  is  a special  case  (U=l)  of  the  one 
in  (36).  As  with  (18),  the  integral  in  (36)  can  be  converted 
into  the  following  integral  in  the  (s,p)  plane: 


p dsdp 


f /_2 


/s^+ (p- ly' 1 ) ^ 


2-m 


bra  a-lj  j 

p dsdp 


f I . I N 2 

l_/s  + (p-  1 Y ' 1 ) 


2-m 


= 2 


Iy' 


a-l  , 

P dp 


ds 


° (/s^+(p-|y'T)^ 


2-u 


rb  , i-a 

t-  I p“~^dp 


ds 


° (/?7(p-1y'  I )^ 


2-u 


= 2|y 


1 , /-a 

a-1  j 

n dn 

0 


ds 


° (/s^+ly' l^(l-n)^ 


2-u 


b/ly' 


a-1  , 
n dn 


ds 


° (/s^+ly'  |^(n-lP, 


2-u 


(C.l) 


where,  above,  we  made  the  change  of  variables  p = |y'ln.  Letting 
also  s = |y'|(l-n)C  and  s = ly'|(n-l)5  in  the  first  and 
second  integrals,  respectively,  we  obtain 


l_/s  +(p-|y'  I ) J 


a-1 


a-1,  „u-i  ra/|y'l(n-l) 
n (n-i)  dn 


We  now  distinguish  two  cases. 

Case  1 ; 0 < u £ y £ 1 . 

With  0 < n < 1 we  have  that 


Similarly,  for  1 < n £ b/jy'j, 


-39- 


a/|y' I (n-i) 


j—  < logl^^ 


rijl  ^ "VI  " [l 


We  note  here  that,  since  b > |y'j  and  because  of  (46),  the 
argument  of  the  first  logarithm  in  each  of  the  last  two 
equations  is  greater  than  one. 


Substituting  (C.3)  and  (C.4)  in  (C.2)  we  obtain 


“"1  j j 

p dsdp 

[/s^+(p-|y'  j ) ^ 


_ < 2|y|“"^-M[  n“-^(i-n)^-Hog 


fk5/|y'  1 p_l 

+ n (n-1)  log 

1 


+ log  1 ^ ^ 9°^  ^(1-n)^  ^dn  + 

I V laj  J Mq 


a-l, 

+ n (n-i)  dn f 

^1 


We  now  assume  that  a+u  > 1.  For  the  first  integral  on  the 
right  (call  it  , we  write 


a+u-2 


h i J n"  . J (l-n)“ 


a+p-2 


dn  = 


(^)a+u-2 


(C.6) 
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I 

f.  I 

f : 


r 


I 

L' 


Case  2 : 1 < u < 2 


With 


y=l+v,  v>0, 


equation  (C.2)  becomes 


p dsdp 


/s2+(p-ly'  1)2^ 


2-m 


= 2ly 


. la+v 


ct-1 ,1  \ V , 

n (l-n)  dn 


0 


a/ly' 1 (i-n) 


(C.IO) 


dl. 


/l+c^. 


1-v 


b/|y' 


n (n-i)  dn 


a/|y’ I (n-l) 


dr. 


i-v 


y < 


< 2ly|“^^ 


a-l , , . V j 

n (i-n)  dn 


a/|y ' 1 (i-n) 
0 


»-l,  ,,v, 

n (n-i)  dn 

1 


a/ly'  I (n-l) 


,1-v 


'3 


it 


2My|“ 


I'o 


a-lj  i 
n dn  + 


J 


b/|y' 

1 


a-l.  2a%“ 

^ = -W- 


(C.ll) 


We  note  that  for  y > 2 the  original  integral  is  well-behaved. 
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Appendix  D;  Proof  of  the  inequalities  in  Equation  (53) . 


For  the  curve  C defined  by  (50)  as 


C = {t  = p+iF(0,p)  = 0 < p < b} 


we  define  the  position  vector  r(p)  by  means  of  (51) 


r(p)  = pp+F(0,p)z  = r(cos  9 p + sin  0 z)  , 


(D.l) 


where  p and  z denote  unit  vectors  in  the  p and  z directions, 
respectively.  We  can  then  write  for  the  tangent  vector 


A+F'^ 


dr  = (p+F'z)dp  = (cos'Jjp+sim})z)  /1+F'  dp  = (cosi|jp+simi)z)  ds  (D.2) 


where  s denotes  arclength  along  C,  and 


cos  = 


1 . F ' 

■ , sin  = 


/l+F'  ^ 


/T+F' 


F ' = 


9F(Q,p) 

9p 


(D.3) 


On  the  other  hand. 


dr  = — dp  = _p+F  F — — cQsi|)  + sin0  sini|j)ds 


/p^+F^ 


/p^+F^/l  + F'^ 


(D.4) 


For 


0 < p < b 


and  b sufficiently  small  we  can  have 
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COS  8 ^ , I sin  6(  £ j 


(D.5) 


since  C is  a Lyapunov  curve  (cf.  Gunter,  1967).  The  same 
statement  can  be  made  about  the  angle  \p  which  is  the  angle 
the  tangent  to  C at  a point  makes  with  the  p-axis  or, 
equivalently,  the  normal  makes  with  the  z-axis.  We  then  have. 


dr  £ (cos0  cosij;  - sin  I 6 1 sin  1 if)  1 ) ds  > 


2 


y'l'  1 


2 2 


ds  = i ds 


(D.6) 


or 


dr  ^ i ds  > ^ dp. 


(D.7) 


To  prove  the  second  of  (53)  we  apply  the  mean  value  theorem 
to  the  function  F: 


= /p^+F(0,p)^  = P /l+(F'  (0,p^)  )^  = 


< — P, 


cos  - ^3 


(D.  8) 


where  0 < p^  < p, 


i. 


A.  Introduction 


One  of  the  most  prominent  problems  of  electromagnetic 
scattering  theory  is  that  of  scattering  of  time-harmonic  waves 
by  a perfect  conductor  which  occupies  a finite  region  of 
space  and  is  bounded  by  a closed  surface  S.  Mathematically, 
this  problem  is  a boundary  value  problem  stated  in  terms  of 
Maxwell's  equations,  radiation  conditions,  boundary  values 
of  the  tangential  component  of  the  electric  field  on  the  surface 
S,  and  continuity  properties  of  the  scattered  fields  in  space 
and  on  the  surface.  One  way  of  studying  this  problem,  i.e. 
answering  questions  on  the  existence,  uniqueness,  and  properties 
of  solutions,  is  to  convert  it  to  a problem  in  integral  equa- 
tions where  the  unknown  function  is  the  linear  current  density 
on  the  surface  S.  This  conversion  can  be  accomplished  provided 
the  surface  and  the  current  density  possess  certain  mathematical 
properties  (cf.  Muller,  1969). 

A problem  very  similar  to  the  above  but  possessing  one 
additional  feature  is  that  of  scattering  of  electromagnetic 
waves  by  an  open,  perfectly  conducting  surface  S.  The  distinctive 
feature  is  the  very  fact  that  the  surface  is  open.  As  pointed 
out  in  Heins  and  Silver  (1955) , it  has  been  known  for  some  time 
that  the  conditions  imposed  on  the  closed  surface  problem 
described  above  are  not  sufficient  in  guaranteeing  the  unique- 
ness of  the  solution  of  the  open  surface  problem,  and  that,  in  order 


to  have  uniqueness,  additional  ones  are  needed.  These  conditions 
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I 


vary  from  author  to  author  (of.  Heins  & Silver,  1955;  Jones 
1964)  but,  ultimately,  all  of  them  have  to  do  with  the 
behavior  of  the  scattered  fields  near  the  edge  of  the  open  sur- 
face and,  in  effect,  demand  that  no  sources  of  electromagnetic 
waves  are  induced  on  the  edge. 

Most  of  the  works  mentioned  in  the  last  two  references 
deal  with  plane  surfaces  (or,  equivalently,  with  apertures  in 
perfectly  conducting  plane  screens) , and  for  the  most  part  aim 
only  at  obtaining  order  relations  for  the  scattered  fields  near 
the  edge  of  the  open  surface,  quite  often  at  the  expense  of 
mathematical  rigor.  In  the  present  work  we  propose  to  formulate 
the  problem  of  scattering  of  electromagnetic  waves  by  an  open, 
perfectly  conducting  surface  as  a boundary  value  problem,  then 
convert  it  to  a problem  in  integral  equations,  and  finally 
prove  that  it  can  have  at  most  one  solution.  In  Section  B we 
present  the  class  of  open  surfaces  that  we  will  consider,  some 
coordinates  systems  and  notation,  as  well  as  a brief  review  of 
the  basic  definitions  and  theorems  of  vector  analysis.  In 
Section  C we  define  the  boundary  value  problem,  i.e.  present  all 
the  conditions  that  the  scattered  electromagnetic  fields  must 
satisfy.  Besides  the  usual  ones  required  of  closed  surfaces, 
as  described  above,  these  conditions  include  statements  on  the 
behavior  of  the  scattered  fields  as  well  as  of  the  induced  current 


I 


density  near  the  edge  of  the  surface.  The  scattered  fields  are 
required  to  satisfy  an  energy  condition  near  the  edge 


I 
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which  requires  that  the  energy  enclosed  in  a finite  region  of 
space  is  finite  and  that  it  vanishes  with  the  volume  of  the 
region.  The  induced  current  density  is  required  to  have  a normal 
component  to  the  edge  which  achieves  a finite  limit  at  the 
edge,  while  its  component  parallel  to  the  edge  is  allowed  to 
grow  beyond  bound  in  the  approach  to  the  edge  but  in  a prescribed 
way.  The  current  density  is  also  required  to  have  first  partial 
derivatives  which,  near  the  edge,  behave  in  the  manner  of  its  com- 
ponent parallel  to  the  edge.  In  Section  D we  assess  the  physical 
implications  of  the  requirements  on  the  scattered  fields  and 
current  density  near  the  edge. 

In  Sections  E and  F we  prove  the  equivalence  theorem, 
i.e.  that  the  boundary  value  problem  of  Section  C is  equivalent 
to  a problem  in  integral  equations,  while  in  Section  G we  prove 
the  uniqueness  theorem,  i.e.  that  the  problem  can  have  at  most 
one  solution.  In  Section  H we  offer  some  concluding  remarks,  and 
in  Appendices  A,  B,  and  C some  detailed  computations. 


Preliminary  considerations 


In  this  section  we  will  introduce  the  scattering  surface, 
two  coordinate  systems  associated  with  it,  some  notation,  as  well 
as  define  the  basic  operations  of  vector  analysis. 

The  surface  under  consideration  is  a perfectly  conducting 
surface  S bounded  by  a curve  C.  We  denote  by  S the  closure 


I 

\ 


i 


} 


I 

I 


! 


-A. 


Of  S(S  = s UC)  and  we  require  that  S is  a regular  open 
surface.  The  definition,  properties,  and  local  description 
of  such  a surface  are  given  in  Part  I,  Section  D,  and  for  this 
reason  we  do  not  repeat  them  here.  In  order  to  reach  the 
objectives  of  this  paper  we  will  have  to  use  this  surface  in 
conjunction  with  the  integral  theorems  of  vector  analysis.  To 
this  end  we  need  to  parametrize  the  curve  C with  respect  to  its 
arclength  and  to  also  introduce  two  coordinate  systems. 

From  an  arbitrary  point  M of  C we  measure  arclength, 
s',  along  C and  we  have  the  following  parametric  representation 
for  the  curve  with  respect  to  an  arbitrary  rectangular  coordinate 
system  xyz: 


x=f(s'),  y=g(s'),  z=h(s'),  0£s'£L,  (1) 

where  L is  the  length  of  C.  VJe  note  that  because  of  the 

definition  of  C in  Part  I,  Section  B,  the  functions  g, 

and  h are  twice  continuously  differentiable  in  s';  moreover, 

~ 2 '■  2 ~ 2 

f +g ' +h ' =1.  The  unit  tangent  vector  to  C is  then  given  by 


t = f(s')  x+g(s')  y+h(s' 


(2) 


The  unit  normal  vector  to  S is  now  chosen  so  that  it  is 
positively  oriented  with  respect  to  t.  This  can  be  done  since 
S is  orientable  by  definition.  With  respect  to  the  coordinates 


of  (l.I) ' the  unit  normal  vector  is  chosen  to  be 

/V  A /N 

-F  x-F  y+z 

n = (3) 

JlAF  ^+F  ^ 

X y 

i.e.  we  agree  that  the  local  z-axis  is  positively  oriented  with 
respect  to  t.  For  points  M of  C,  we  call  h of  the 

unit  normal  vector  on  C.  If  we  also  define 

T = n X t (4) 

then  the  triple  (n,t,x)  is  a positive  triple  of  orthonormal 

vectors  at  each  point  of  C.  We  note  that  the  vectors  t and 

T lie  on  the  plane  tangent  to  S at  the  point  of  C under 
consideration . 

The  remaining  coordinate  system  is  a polar  system  ( p ',({)' ) 
p'  ^ 0,  0 £(})'£  2-rr.  It  is  erected  on  the  plane  of  the  vectors 

A A. 

T and  n with  the  pole  at  the  origin  of  this  plane,  and  the 
angle  <1> ' measured  from  x to  ix.  If  r-Q  is  the  position 

vector  to  a point  of  C,  and  if  p'  is  the  position  vector  on 

the  x-n  plane,  then  the  equation 

^ A A A 

r = xx+yy+zz  = r^+p’  (5) 


(1)  The  Roman  numeral  one  (I)  refers  to  Part  I of  this  work. 
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defines  a transformation  of  points  (s',p',(p')  to  points 
(x,y,z).  The  Jacobian  of  this  transformation  is 


3p'  9(j) 


dr  _ , f. 


The  continuity  of  the  second  derivatives  of  f,  g,  and  h 
guarantees  the  boundedness  of  | dt/ds ' | which  in  turn  guarantees 
that,  for  p'  sufficiently  small,  the  Jacobian  is  positive 
and  the  transformation  one-to-one. 


At  this  point  we  have  completed  the  discussion  of  the 
coordinates  systems  that  will  be  employed  below.  Before  moving 
to  the  subject  of  vector  analysis  we  introduce  some  notation  that 
will  be  used  in  the  following  section.  We  denote  by  S_^  the  side 
of  S facing  in  the  direction  of  the  normal,  and  by  S_  the 
other  side.  If  f is  a function  defined  in  a region  containing 
S,  we  denote  by  f^(f_)  its  limit  as  its  argument  approaches 
S from  S^(S_) . Finally,  we  denote  by  the  closed  curve 

which  is  the  intersection  of  S with  the  cylinder 


0 £ s'  £ L,  p'  = G > 0,  0 £ ({)  ' £ 2n, 


where  we  take  g small  enough  so  that  there  is  an  one-to-one 
correspondence  between  points  of  and  points  of  C.  Due  to  the 

smoothness  of  S,  possesses  a tangent  vector  t which  we 

take  positively  oriented  with  respect  to  the  normal.  If 
is  a point  of  such  that  M g C,  as  g -►  0,  then 

t(M^)  ->  t(M)  . 


We  conclude  this  section  by  explaining  the  sense  in  which 
the  basic  operations  of  vector  analysis  are  to  be  understood. 
The  ideas  that  follow  have  been  taken  from  Muller  (1969)  who 
presents  them  in  his  book  in  great  detail  so  that  we  will 
present  here  only  what  is  absolutely  necessary  for  our  work. 

To  do  so  we  assume  familiarity  with  the  definition  of  a regular 
surface  and  a closed  regular  surface  (Kellogg,  1953;  Muller, 

3 

1969) . A regular  region  will  be  a compact  point  set  in  R 
bounded  by  a closed  regular  surface. 

Definition  1:  A sequence  of  regular  regions  is  said  to 

converge  to  a poxnt  r^  if  for  every  e > 0 there  exists  a 
number  N(e)  such  that  all  with  v ^ N(e)  are  entirely 

within  the  region  | r-rQ ] £ e. 

>•  ► 

Definition  2:  Let  v(r)  be  continuous  in  the  neighborhood  of 

the  point  r^.  Denote  by  | | G^ | | the  volume  of  the  regular 
region  G^  and  by  = 9G^  , its  boundary.  If,  for  each 

sequence  G^  converging  to  , the  limit 


1 im 

G ->r 
V 


0 


- - 
n • vdF 
F 

V 


exists  and  is  unique,  then  we  set 


V • V = lim 
G..--r, 


n • vdF , 


where  n is  the  unit  normal  to 


divergence  of  v at  r 


0 


We  call  V*v  the 


Definition  3:  With  the  notation  and  conditions  of  Definition  2, 

if  the  limit 


lim 

G -^r 
V 


0 


nxv  dF 
F 

V 


exists  and  is  unique  for  each  sequence  converqinq  to 

then  we  set 


r 


0' 


V X V = lim 


nxvdF 


G -^r.l  |G 
V O'  I \ 


We  call  Vxv  the  curl  of  v at  r^ . 

These  definitions  differ  from  the  usual  ones  in  terms  of 
the  operator 


^ ^ A rv 


in  that  they  do  not  require  v to  have  first  partial  derivatives 
can  be  shown,  however,  that  if  v has  continuous 
first  partials  in  the  neighborhood  of  r^ , then  the  two 
definitions  are  equivalent  (Muller,  1969).  From  the  definitions 
above  we  have  the  following  two  theorems  (Muller,  1969). 


Theorem  1:  Let  v be  continuous  in  the  regular  region  G. 

Let  V'v  be  continuous  in  each  subregion  lying  entirely 
If  the  integral 

[ V • V dV 


i 


! 


I 

1 

i 


in  G. 


where  F = 9G. 


Theorem  2:  Let  S be  a surface  bounded  by  the  curve  C so 

i. 

— ->•  , 
that  S = SVJC  is  a regular  open  surface.  Let  v and  Vxv  M 

h 

I j 

be  continuous  in  a region  that  contains  C.  Then  ' 


n*Vxv  dS  = 
S 


These  two  theorems  are  the  well-known  divergence  and 
Stokes'  theorems,  respectively.  It  must  be  noted  that  they 
do  not  require  differentiability  of  the  vector  v.  Besides  these 
theorems  we  will  need  the  surface  divergence  theorem  and  for 
this  reason  we  introduce  the  following  definitions 


t • V ds ' . 

■*C 


Definition  4 : A sequence  of  regular  open  surfaces 

U is  said  to  converge  to  a point  rg  if  for  every 

e > 0 there  exists  a number  N(s)  such  that  all  S with 

V 

• • 1 I 

V ^ N(e)  are  entirely  within  the  region  | r-rg | ^ e. 

Definition  5:  Let  v(r)  be  continuous  in  the  neighborhood  of 

the  point  r^.  Denote  by  | | | | the  surface  area  of  the 
regular  open  surface  S^.  If,  for  each  sequence  converging 

to  r,, , the  limit 


-11- 


exists  and  is  unique,  then  we  set 


VU  = lim_^ 

G 

V 0 


n UdF , 


and  we  call  VU  the  gradient  of  U at  r. 


The  open  surface  problem. 


In  this  section  we  present  in  mathematical  form  the 
problem  of  scattering  of  electromagnetic  waves  by  a perfectly 
conducting  regular  open  surface  S = S C . Before  doing  so 
we  introduce  a definition  that  we  will  need  below.  In  Part  I, 
Section  C,  we  defined  a class  of  functions  that  we  called 
regular  density  functions.  We  now  restrict  this  class  to  the 
following. 

Definition  7.  A real-valued  function  g defined  over  a regular 
open  surface  S is  said  to  be  an  H-regular  density  function  for 
the  surface  if  it  satisfies  Eq.  (1.14)  and  the  following  condition; 


In  every  closed  and  connected  subset  S'  of  S,  bounded 
away  from  C,  the  function  g is  Holder-continuous,  i.e.  there 
exist  real  numbers  B and  3 such  that 


g (Mj^) -g  (M^)  1 < b|Mj^-M21®,  0 < 6 £ 1, 


where  B does  not  depend  on  the  position  of  and  M2  in 

S',  but  may  depend  on  the  proximity  of  S'  to  the  boundary 
C of  S. 
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l 

I 


i 

i 

I 


i 


i 

\ 


The  open  surface  problem  can  now  be  stated  as  follows; 

We  wish  to  find  vector-functions  E and  H which  are  defined 
3 3 — 

in  R -C,  are  continuous  in  R -S,  also  continuous  to  the 
surface  S from  and  S_,  and  satisfy  the  following 

conditions 

(i)  Maxwell's  equations 


VxE(R)  = ikZH(R),  Vxh(R)  = -ikYE(R),  R e (R^-S) 


(8) 


where  Z and  Y are  the  free-space  impedance  and  admittance, 
respectively,  so  that 


z . v-i  = , 


(9) 


with  Pq  the  permeability  and  the  permittivity  of  free 

space.  Equations  (8)  are  also  required  to  hold  in  the  approach 


to  S from  and  S 


(ii)  The  Silver-Muller-Wilcox  radiation  conditions 


/N  , 

— ► 

fll 

'i' 

+ ZRxH  = o 

, H - YRxE  = o 

R 

R 


(10) 


uniformly  in  the  non-radial  directions.  Here,  R = RR,  with  R 
the  unit  vector  in  the  direction  of  R , 


(iii)  The  edge  conditions 


I 


(a)  For  M e C and  with  respect  to  the  polar  coordinate 
system  of  Section  B 


-13- 


a > 0 , 
(11) 


2it 

E*E*d<j) ' 

0 


0 (p'^ 


f 27T 

H-H*d(f) ' 

■I  0 


where  the  order  is  uniform  with  respect  to  M e C, 
and  where  the  asterisk  denotes  the  complex  conjugate 
of  the  function  it  is  attached  to. 

(b)  On  a family  of  curves  ->  C as  described  in 

Section  B,  the  function 

= t(M^)  • (H_  (M^) -H^  (M^)  ) , (12. a) 

A 

where  M e C and  t(M  ) the  unit  tangent  vector 
e e e 

to  at  , has  a finite  limit  as  e -*•  0 , 

and  the  limit  function  i|^(M),  M e C,  defined  by 


ij;(M)  = lim  4)  (M  ),  (12. b) 

e^O  ^ 

with  ->■  M as  e -*•  0 , is  continuous  at  every  point 

M of  C.  Moreover,  the  function  4*  defined  by 


4^  (M^) 


i[)(M^)  , 0 < e £ e ' 

' b(M) , e = 0 


(12. c) 


is  an  Holder-continuous  function  of  c on  0 £ c £ ^ ' 
where  e'  is  small  enough  so  that  there  is  an  one-to- 
one  correspondence  between  points  of  and  points 


of  C. 


(iv)  The  boundary  condition 


nxE 


= nxE 


= -nxE 


(13) 


where  is  the  electric  field  of  a specified  plane  wave  or  a 

dipole  source  or  a combination  thereof  located  off  S.  (The 
condition  on  the  normal  component  of  the  magnetic  field,  i.e. 


->■ 

-y  , 

n • H 

= n*H 

= -n*H 
S_ 

need  not  be  specified  since  it  can  be  derived  from  (13)  by 
taking  the  divergence  of  it.) 

(v)  The  density  conditions:  that  the  function 

k(r)  = n(r)  x (H_ (r) (r) ) , r e S,  has  partial  derivatives  of 
the  first  order  on  S,  and  that  the  real  and  imaginary  parts  of 
the  components  of  k as  well  as  those  of  its  first  partials  are 
H-regular  density  functions  for  S (Def.  7). 

■ --  With  respect  to  this  statement  of  the  open  surface  problem 
we  would  like  to  make  the  following  comments'.  the  continuity 

of  £ and  H in  space  and  in  the  approach  to  S from  either 
side  together  with  condition  (i)  imply  that  Vx£  and  Vxh  are 
continuous  in  space  and  in  the  approach  to  S from  either  side. 
Second,  condition  (iii.b)  is  needed  only  in  proving  the  uniqueness 
of  the  solution  of  the  open  surface  problem  but  not  for  converting 
it  into  a problem  in  integral  equations.  The  same  is  true  about 
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the  existence  of  the  first  partials  of  K in  condition  (v) . 
Indeed,  for  the  conversion  part  we  will  only  need  that  K and 
Vq*K  (see  Def,  5)  have  real  and  imaginary  parts  which  are 
H-regular  density  functions  for  S.  As  condition  (v)  stands 
Vq*K  exists  and  is  equal  to  the  one  given  in  terms  of  the 

differential  operator  V since,  by  definition  7,  the  first 

• • • 

partials  of  K are  Holder-continuous  on  S. 

The  problem  described  above  differs  from  the  corresponding 
closed  surface  problem  mainly  in  the  addition  of  conditions  (iii) 
Before  undertaking  the  main  task  of  this  paper,  i,e.  to  convert 
this  problem  to  a problem  in  integral  equations  and  to  also  show 
that  it  has  at  most  oi.  solution,  we  would  like  to  first  discuss 
the  edge  conditions  and  their  physical  implications. 


D.  The  edge  conditions. 


It  is  well-known  that  in  order  to  guarantee  the  uniqueness 
of  solution  of  the  open  surface  problem  we  need  one  more  condition 
than  what  is  required  for  the  closed  surface  problem  (Jones,  1964; 
Sommerfeld,  1964. b).  This  condition  is  associated  with  the  elec- 
tromagnetic  energy  enclosed  in  a finite  region  of  space.  If  e 
and  h are  electromagnetic  fields  satisfying  Maxwell's  equations 


vxe  - -Ko  t , VxR  - If- 


(15) 


where  Eq  and 


stand  for  the  permittivity  and  permeability 


I 
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of  free  space,  respectively,  we  require  that  the  stored  electric 
and  magnetic  energy  in  a finite  region  of  space  devoid  of 
sources  be  finite.  For  a region  V of  space  these  energies 
are  (Sommerfeld,  1964. a) 


W = -j  I e-e  dV,  w = ~ [ 
e 2 m 2 


h'h  dV. 


If  e = Re (Ee  ^) , h = Re (He  ) , with  E and  H depending 
on  the  space  variables  only,  then  the  time-average  electric  and 
magnetic  energies  are  given  by 


W = E-E  dV,  w = ~ H-H  dV, 

e 4 m 4 

J V ■'  V 


respectively.  The  bar  denotes  time-averaging  over  one  period, 
while  the  star  denotes  the  complex  conjugate  of  the  vector  it 
is  attached  to  (Jones,  1964;  Muller,  1969;  Stratton  , 1941). 

In  the  present  case  the  volume  would  be  a tubular  one 
surrounding  the  edge  C of  S.  With  respect  to  the  coordinate 
system  introduced  in  Section  B,  this  volume  (call  it  V^)  is 
defined  by 


0 _<  s'  < L,  0 £ p'  £ e,  0 £ 4>'  £ 2tt 


The  finiteness  of  energy  requirement  can  then  be  written  as 


lim  E-E*dV  = 0, 
€;--0 


lim  H-H  dV  = 0 
c-^0 


Instead  of  this  requirement,  we  impose  the  edge ■ conditions 


•2it 


E*E  d4> 


0 (p 


,2 (a-1) 


f2TT 


i-H*d(j)’  = ) , a > 0 


where  the  order  conditions  are  taken  to  be  uniform  with  respect 
to  the  point  M of  C.  Clearly,  if  a pair  of  functions  E and 
H satisfy  the  edge  conditions,  then  they  also  satisfy  the 
energy  conditions. 

Although  conditions  (20)  are  necessary  for  proving  that  the 
open  surface  problem  has  at  most  one  solution,  they  are  not 
sufficient,  at  least  not  for  the  manner  in  which  we  approach  the 
question.  Condition  (iii.b)  of  Section  C is  also  needed  and  this 
we  will  demonstrate  later  on.  In  the  meanwhile,  we  will  show  that 

the  introduction  of  such  a condition  makes  sense  physically.  Let 

->■  -► 
f(R)  be  a complex-valued  function  of  R defined  and  with 

continuous  first  partials  in  a bounded  domain  containing  S in 

its  interior,  and  let  E and  H be  a solution  pair  for  the  open 

surface  problem.  Let  also  C^  represent  the  closed  contour  on 

S which  is  the  intersection  of  the  surface  S with  the  cylinder 

in  (7).  By  Stokes'  theorem  we  have  that 


f (R)t 


(H_-H^)ds 


f _ _ 


f (R)  t'fl  ds  ' - 


-C 


f (R)fH^ds' 


r - 


n*Vx (fH)dS 
e 


[7f  xii+fVxH]dS 


[Vf • (pxH) +ikYfp ' -E]dS, 


(21) 
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where  -C^  is  transversed  in  the  opposite  direction  of  C^, 


s'  denotes  arclength,  denotes  the  cylinder  in  (7),  and 

n the  unit  normal  to  the  cylinder  pointing  toward  the  interior. 


We  can  now  show  that  the  last  integral  in  (21)  vanishes  in 
the  limit  as  c 0 . Since  E and  H are  solutions  of  the 
open  surface  problem  they  satisfy  (20).  Using  the  Cauchy- 
Schwarz  inequality,  we  have  that 


fp' *EdS  ^ < 1 f I I p' -E|dS  ^ 


L .2tt  ^ ^ ^ ^2  ) 2 

|f  I Ip-  -Eje  l-ep'  • ^ d4>'ds'  < 

0 J 0 J } 


[ rh  r2T\  , . 2 .L  r2TT 

< e If  r 1-ep'  • d(()'ds' 


p ' • E I ^d(() ' ds  ' [■  < 


where  is  a constant.  Similarly, 


Vf  • (pxft)dS  £ 


where  C2  is  a constant.  Since  a > 0,  both  of  these  expressions 
go  to  zero  with  e.  We  then  have  that 


lim  f(R)t  • (H  -H  )ds'  = 0. 
e-^0  ■'C 

£ 
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But  the  integrand  in  this  expression  is  continuous  to  the  boundary 


by  the  edge  condition  (iii.b)  and  the  assumption  on  f.  This 
allows  us  to  take  the  limit  before  integrating,  so  that 


f (R)  t • (H_-H_|_)  ds  ' = 0 


In  order  to  interpret  these  last  two  equations 

-► 

physically  we  introduce  the  induced  linear  current  density,  K, 
on  the  surface  S 


K(r)  = n(r)  x (H_  (r ) (r ) ) , r e S 


a vector  tangent  to  the  surface  S.  Since, 


t • (H  -H^)  = (Txn)  • (H  -fp)  = T-n  x (h  -5  ) = x-K  (27) 

“ T*  “ *T*  “ T 


we  have  from  (24)  and  (25) 


lim  1 f(R)x'Kds'  = 0,  ] f(R)T’Kds'  = 0. 

£'>•0  Jc  -IC 


With  f(R)  = 1,  and  since  t-K  is  normal  to  C,  the  second  of 
(28)  says  that  the  current  (in  amperes)  entering  the  surface  must 
be  equal  to  the  current  leaving  it.  The  first  of  (28)  says 
basically  the  same  thing.  Physically,  however,  no  current  is 
entering  or  leaving  the  surface.  In  fact,  we  should  have  that 
x*K  = 0 on  C.  This  information  is  contained  in  the  second  of 
(28)  but  not  in  the  first.  To  prove  it  we  use  the  following 
theorem  (Sobolev,  1964,  p.  113). 
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Theorem  4:  If  a function 

and  if  for  any  function  f 


equality 


fgdv  = 0 holds, 

0. 


g is  summable  in  an  open  set  Q, 
which  is  continuous  on  Q the 

then  g must  satisfy  the  condition 


I g I dv  = 0,  and  consequently  g is  equal  to  zero  almost 
Q 


everywhere  in  T2 . (We  remark  here  that  a much  stronger  version 
of  this  theorem  can  be  found  in  Smirnov  (1964,  p.  145)). 

Identifying  g with  t*K  in  the  second  of  (28),  we  have 

/V  ^ 

that  T’K  is  equal  to  zero  almost  everywhere  on  C.  Since  by 

^ A. 

(27)  and  (iii.b)  x*K  is  continuous  on  C,  then  t-K  = 0 every- 
where on  C.  Thus  the  edge  conditions  of  the  open  surface  problem 
imply  that 

fit  = 0 on  C.  (29) 

The  second  of  (28)  or,  equivalently,  (25)  provides  us  with 
another  physical  statement.  Using  the  surface  divergence 
theorem  (theorem  3) , we  have  that 


or. 


0 


f (R)T- Kds'  = -f  V • (fK)dS, 
C J S 


! 

i 

- s 


V 


0 


(fK)dS 


0. 


(30) 


With  the  accepted  definition  for  the  electric  charge  density, 
a,  induced  on  S,  i.e. 


a 


£gn • (E_-E^) , 


(31) 
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and  the  fact  that 

/\ 

V„-K  = -n*V-x(H  -H^)  = ikYn* (E  -E^)  = iwo,  (32) 

we  have  from  (30) , with  f = 1,  that 

ads  = 0,  (33) 

S 

which  says  that  the  electric  charge  induced  on  S is  equal  to 
zero . 

E . The  integro-dif ferential  equation  and  integral 
representations . 

In  this  section  we  convert  the  open  surface  problem  of 
Section  C into  a problem  in  integral  equations. 

>•  > 

Theorem  5 : If  the  pair  of  vector-functions  {E,H)  is  a solution 

-y  •> 

of  the  open  surface  problem,  then  E and  H are  given  by 


E(R')  = 


Z 

ik 


V •K(R)VG(R|R' )+ikZG(R|R' )K(R) 


dS, 


R' 


(R^-S) 


(34) 


H(R’)  = 


K(R) xVG(R[R' )dS,  R'  e 


(R^-S) 


(35) 


^ik|R-R' I 
4tt|  r-R'  I 


where , 


G(r| R' ) 


R 5^  R'  , 


and  where  the  linear  current  density  K is  defined  in  (26)  and 
satisfies  the  integro-dif ferential  equation 


-n(r)xE-^(r)  = 


VQ-Kn(?)  xVG(R|r)+ikZG(R|?)n(?)  xK(r) 


dS , re 
(37) 


Moreover,  on  a family  of  curves  -»■  C,  the  function 

t(M^)*K(M^),  e , is  identical  to  the  function  \[j  defined 

in  ( 12 . a) , and 


limT(M)‘K(M)=0.  (38) 

e->0  ^ ^ 


Proof ; To  show  that  E is  given  by  (34)  we  surround  the  surface 

5 by  a closed  surface  Sq  which  we  construct  as  follows:  Let 

6 and  e be  two  positive  real  numbers  such  that  0 < 6 < e. 

If  Rg  is  the  position  vector  describing  S,  we  let 


S 


+5 


6n  , 


(39) 


Moreover,  we  let 


I 

e 


e . 


(40) 


The  surface  Sq  is  the  closed  surface  formed  by  these  three 
intersecting  surfaces  as  shown  in  Figure  1. 


We  apply  now  Green 
bounded  by  the  surface 


second  identity  in 
the  surface 


the  volume  Vq 
of  a sphere  whose 


radius  R eventually  recedes  to  infinity,  and  the  surface  S 


With 


v/here  G is  given  by  (36),  and  I is  the  identity  dyadic,  we 


have  that 


where  the  unit  normal  n 


to  S„  points  away  from  V 


both  H and  F 


the  volume 


Kleinman,  1971) 


(2)  Double  bars  over  a letter  denote  a dyadic.  All  dyadic  and 
vector  identities  used  in  this  work  can  be  found  in  Van  Bladel 
(1964)  . 


(3)  This  identity  can  be  obtained  by  writing  the  dyadic  in  (41) 
in  rectangular  components  and  applying  theorem  1.  The  vector  to  be 
used  in  this  theorem  is  of  the  form  v = u x w.  From  Muller  (1969) 

-V  -»•  - ->• 

we  have  that  if  u,  w,  V^u,  and  Vxw  are  continuous  in  , 


These  conditions  are  satisfied  in  (42) 
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r->0  •>  S 


lim  [HxVxr+(VxH) xp]dS  = -ikV'^H(R') 


By  these  last  two  equations  and  (8)  we  have  that 


E(R')  = -y  [ n„ • [H(R) xVxr (r(r' )+yxH (R) xr (R  R' ) ]dS.  (44) 

J q +q 


Vxr  = ik[VVG+k  GI]  , 


HxWG  = (VxH)  VG-Vx  (hVG) 


then,  denoting  the  integral  over  by  I,  we  have  that 

K 


1 = i |-ikR-EVG+k^GZRxlS-ik  [R-VGE-RE*  VG]  jdS 


Employing  polar  coordinates  with  center  at  R'  and  writing 


- ^ e^’^^R(ZRxH+E)dn  + e^’^^SdQ. 


The  first  of  these  integrals  vanishes  in  the  limit  as  R 


We  next  show  that  the  integral  over  vanishes  in  the 

K 

limit  as  R ->  «.  Since, 


e^^^  fik  ll  ''  7 

= - — -R  - = R^dfi 

R 


we  obtain 
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because  of  (10),  while  the  second  vanishes  because  of  a lenuna 


of  Wilcox  (1956)  which  in  effect  states  that  Maxwell's  equation 


together  with  the  radiation  conditions  imply  that 


I I 2 

E d(2 


R 


Thus  (44)  becomes 


E(R')  = I n- [H(R)xVxr(R|R')+VxH(R)xr(R|R') ]dS, 

K bQ 


=-►  ->.  — - 


which  by  (45)  and  (46)  can  be  written  as 


S(R')  = n„*  [Hx(VVG+)c^GI)-il<YExVx(GI)  ]dS. 

K J U 

^0 


In  order  to  cast  this  expression  in  the  usual  Stratton-Chu 


form  we  use  the  identities 


n^ • (Hxl)  = nxH 


n_-  [ExVx(GI) ] = n- • (Ex(VGxl) ] = (n„ xE)  • ( VGxi ) = (n„xE) xVG  (54) 


nQ-(HxVVG)  = (nQXH).VVG  = -(nQXH)-V'VG  = -V'(nQH*VG)  = 


= V'(nQ-VGxH)  = V [nQ-Vx  (GH)  ]-V  (GnQ-VxH)  = 


= V (Hq -Vx  (GR)  ] -iJcYnQ -EVG  . 


(55) 
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Eq.  (52)  then  becomes 


E(R')  = 


[nQ-EVG+ikZGnQxH+(nQxE)  xVG]dS  + nQ-Vx(GH)dS. 

^0  ^0 


The  last  integral  in  this  expression  vanishes  by  virtue  of 
Stokes'  theorem.  Letting  now  6 ->■  0 , and  using  the  boundary 
condition  (13),  we  have  that 


E(R')  = 


[nQ-EVG+ikZGnQxH+(nQXE)'xVGldS  + 


+ [n-(E  -E_^)  VG+ikZGnx  (H  -H^))dS, 

Js  - + - + 

e 


where  S is  that  portion  of  S not  contained  in  E . The 
first  of  these  integrals  can  be  shown  to  vanish  as  e -*■  0 by 
means  of  the  edge  conditions  (11) . Since  R'  / , G and  VG 

are  continuous  functions  of  R and,  hence,  bounded.  Using  the 
Cauchy-Schwarz  inequality  as  for  (22) , we  have  that 


A ^ 

Uq  -EVGdS  < )C^, 


L - 


qXE)  xVGdS  £ (R'  ) e 


GnQxHdS  £C2(R')e 
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so  that  the  first  integral  in  (57)  is  indeed  zero  in  the  limit 
as  £ ->■  0 . The  second  integral  in  (57)  can  be  shown  to  exist 
as  an  improper  integral  by  condition  (v)  of  the  open  surface 
problem.  For  K = n‘(H_-II^)  is  integrable  since  its  components 
have  real  and  imaginary  parts  that  are  H-regular  density  functions 

A,  - 

for  S.  From  (32),  n*(E_-E^)  = (ikY)  and,  because  of  the 

conditions  (v)  on  K,  exists  and  its  real  and  imaginary 

parts  are  H-regular  density  functions  for  S,  hence  integrable. 
Letting 

’ r 

( • ) dS  = lim  ( • ) dS  , 

■’  S e-^0  J S 

e 

we  have  that 


Vq-K(R) VG(R|R')+ikZG(R|R')K(R)  dS , R'  £ (R^-S) . 


This  completes  the  proof  of  (34).  The  proof  of  (35)  follows 
similar  steps  and  for  this  reason  we  omit  it. 


Equation  (38)  and  the  identification  of  t*K  with  \p  were 
proven  in  Section  D so  that  it  only  remains  to  prove  (37) . This 
equation  can  be  obtained  either  from  (56)  or  (62) . In  (56)  the 
surface  Sg  is  a closed  surface,  a fact  which  enables  us  to 
use  well-established  results  regarding  the  behavior  of  the 
integral  as  R'  approaches  or  S_^  (cf.  Muller,  1969, 

Ch.  IV).  In  (62)  the  surface  is  open  but  this  equation  involves 
only  two  terms  and  for  this  reason  it  is  the  equation  we  will  use 


r 
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to  obtain  (37).  The  theorems  necessary  for  the  derivation  are 
summarized  in  Appendix  A.  From  theorem  A. 3 we  have  that  the 
second  part  of  the  integral  in  (62)  is  continuous  at  every 
point  r E S so  that 


K(R)G(R|  R'  )dS  = 


K(R)G(R|R')dS  = 


K(R)G(Rlr)dS  (63) 


where  by  S^(S_)  we  mean  the  limit  of  the  integral  as  R' 
approaches  r from  the  side  S^(S_) . From  theorem  A. 4 we  have 
that  the  first  part  of  the  integral  in  (62)  is  continuous  at 
every  point  r of  S in  the  approach  to  it  from  either  side 
or  the  side  S_  and  that 


(S_) 


Vq-K(R)VG(R|R' )dS  = 


+ I Vq  *K  (r)  n (r)  + 


V--K(r)n  ^ G(R|r)dS  + 
U on 


Vq-K(?)  VQG(Rj'r)dS  + 


[VQ-K(R)-VQ-K(r) ]VG(r| r)dS,  (64) 


where  the  minus  sign  corresponds  to  the  approach  from  and 

the  plus  from  S . Combining  the  last  three  equations  we  get 


A 

E^(r)  = + VQ-K(r)  + ~ • K (R)  VG  (R  | r ) tikZG  (R  | r ) K (R) 


dS,  reS, 
(65) 


J 
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where  the  first  part  of  the  integral  must  be  interpreted  in 
the  sense  of  (64).  This  last  equation  together  with  the 
boundary  condition  (13)  yield  (37) , and  the  proof  of  Theorem  5 
is  complete. 


F.  The  converse  of  Theorem  5. 


Before  we  state  and  prove  the  converse  of  Theorem  5 we 
introduce  the  following  definition 

Definition  8.  Let  K be  a vector-function  defined  on  S.  Then 
K IS  said  to  be  a surface  field  on  S if  n-K  = 0 everywhere 
on  S . 


Theorem  6.  Let  K be  a surface  field  defined  on  S and  satisfying 
the  following  conditions 

-► 

(i)  K has  partial  derivatives  of  the  first  order  on  S,  and 

— W 

the  real  and  imaginary  parts  of  the  components  of  K as  well 
as  those  of  its  first  partials  are  H-regular  density  functions 


for  S (Def . 7) . 


(ii)  On  a family  of  curves  ->■  C as  described  in  Section  B, 

the  function  li^(M^)  = x(M^)’K(M^)  has  a limit  as  e 0 , 


lim  iJ^(M  ) = 0 
c->-0 


;66 .a) 


Moreover,  the  function  4*  defined  by 


ip  (M^ ) , 0 < e £ e ' 

0 , c = 0 


(66. b) 
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is  an  Holder-continuous  function  of  e on  0 £ ^ 

where  e ' is  small  enough  so  that  there  is  an  one-to- 

one  correspondence  between  points  of  and  points  of  C. 


(iii)  For  all  r e S,  K satisfies  the  equation 


-n(r)xE^(?)  = VQ-K(R)n(r)  xVG(Rir)+ikZG(R|r)n(r)  xk(R)  dS, 

S 


where  is  the  electric  field  introduced  in  condition 


(iv)  of  the  open  surface  problem. 


Then  the  functions  E and  H defined  by 


E(R')  = •K(R)VG(r1r' )+ikZG(R|R' )K(R)  dS,  R'  e (R^-S) 

Js  ^ J 


H(R')  = K (R) xVG(R| R' ) dS,  R'  e (R  -S)  (69) 

J S 


where  G is  defined  in  (36),  are  solutions  of  the  open  surface 
problem. 

To  facilitate  the  proof  of  this  theorem  we  first  introduce 
the  following  lemma. 


Lemma  1:  If  the  hypotheses  of  Theorem  6 hold,  then 


[G(R|R')V^*K(R)+K(R) •VqG(R|R’) ]dS  = 0,  R'  c (R^-C) . (70) 


Proof:  For  R'  e (R'^-S)  we  have,  by  means  of  the  surface  diver- 


gence theorem  (Theorem  3)  and  (66),  that 


i 
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0 = G{R1r'  ) x-K(R)ds'  = 

Jc 


- V • [g(r|r')k(r) ]dS 
Js 


” 3 — 

Expanding  the  surface  integral  we  obtain  (70)  for  R £ (R  -S) . 

From  Theorem  A. 3 we  have  that  the  first  part  of  the  integral  in 

(70)  is  continuous  at  every  point  reS.  The  same  is  true  for 

the  second  part  and  to  prove  this  we  erect  a rectangular  coordinate 

—V 

system  xyz  with  origin  at  r and  the  z-axis  in  the  direction 
of  the  normal.  We  can  then  write 


K(R)  -V  G(R|R' )dS  = K (R)~  + K (R)—  + K (R)|f  dS  = 
g 0 ' X dx  y 3y  z dz 


= X-  K (R)VGdS+y  K (R)VGdS  + z-l  K (R)VGdS. 
Jq-X  J c;  y 


= ikz[  K(R) xVG(R|r' )ds  = ikZH(R'). 
■’  S 


V'xH(R')  = V X [K(R) xVG(RlR' ) ]dS  = 
^ S 


[K(R)V  G(R|R')+K(R) •V'VG(rIr') ]dS. 


Applying  the  results  of  Theorem  A. 4,  with  n(r)  = z and 
K^(r)  =0,  to  this  equation  we  get  the  continuity  property  at 


r e S . 


We  now  proceed  to  prove  Theorem  5.  We  first  show  that  (68) 

-I-  3 — 

and  (69)  satisfy  Maxwell's  equations  (8)  at  all  points  R'  e (R  -S) 


V'xE(R')  = ikZ  V'x[G(r|R' )K(R) ]dS  = 
S 
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By  a standard  identity  and  (70) 


K(R) -V'ycCRlR' )dS  = V 


K(R)  •VG(r1r' )dS  = -yj  G(RjR’ )VQ-K(R)dS, 
S J S ^ 

(74) 


so  that  (73)  becomes 


V'xH(R')  = -f  [V  •K(R)VG(r1r' )-k^G(R|R' )K(R) ]dS  = ikYE(R'). 

J S 

(75) 

— V ^ 

We  note  that  in  (74)  we  also  used  the  fact  that,  since  K is 
a surface  field,  K*VG  = K*VqG. 

The  functions  E and  H are  clearly  infinitely  differen- 
tiable at  every  point  R'  e (R^-S) , From  (68)  and  Theorems 
A. 3 and  A. 4,  E is  continuous  to  the  surface  S from  and 

S_.  From  (69)  and  Theorem  A. 5 the  same  is  true  of  H.  Moreover, 
from  (72)  and  (75)  and  Theorem  A. 3 and  A. 5,  we  see  that  Maxwell's 
equations  hold  on  S in  the  approach  to  it  from  and  S_. 

We  next  show  the  radiation  conditions  (10)  are  satisfied. 
Asymptotically , 


G(RrR')  = - 


where 


ikR' 


R' 


+ 0 


1 1 

• ' V- 


Uk<t)R’)  + 0 


, R'->oo, 

(76) 


-ikR' *R 
4it 


(77) 


Equations  (76)  are  uniformly  valid  for  all  R c S and  all 
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R' . From  (68)  we  have  that 


E(R')  = [-4)R' VQ-K+ikct)K]dS  + 0 


while  from  (69) 


R'  xii  (R-  ) = - 


i)c(t)R' X (KxR' ) dS  + 0 

. /-I  y-Y  1 ^ 


ik<{)  [K-R'  (R‘  -K)  ]dS  + 0 , 

^ Js  R'^ 


so  that 


Z p 

E (R' ) +ZR' xH (R' ) = R'  — 


-K+ikifiR' -K]dS  + 0 -~ 
^ Js  R' 


We  now  note  that 


7 f— 

V<j)  = - ^V[e  = -ik4>V(R'-R)  = -ik<|)R' , 


so  that  (80)  becomes 


E(R' ) +ZR' xh(R' ) = -R' 


7 / -V  I 

' [(})V.-K+K-V(})]dS  + 0 

^ Js  ^ R' 


-R'  — V -(<))K)dS  4-  0|-^ 

R'  ^ iR’ 


' I " - 1 II 

())T*Kds  + 0 — ^ = 0 — , 

^ J C R'  R’ 


where,  above,  v/e  used  the  surface  divergence  theorem,  eq.  (66), 

« 

and  the  fact  that  K*V$  = K-Vq4>.  Similarly, 


H(R')-YR'xe(R')  = - ik(t)KxR'dS  - 

^ J c: 


f - ( I ] f 1 1 

ik(J)R'xKdS  + 0 = 0 

^ Js  R'2  r.2 


We  note  that  the  first  integral  in  (82)  is  the  first  term  in  an 
asym,ptotic  expansion  of  (70)  . Indeed,  the  radiation  conditions 
can  be  obtained  by  first  using  (70)  and  then  expanding 
asymptotically.  This  we  do  in  Appendix  B. 

We  next  prove  that  the  edge  conditions  are  satisfied.  To 
this  end  we  let 


H(r|r' ) = G(R|R' ) + 


4tt  R-R' 


“ - kl"''  + |R-R'|  + . 


We  can  then  write  in  place  of  (68) 


Vg-K(R)VH(R|R')+ikZH(R|R')K(R)  dS,  R'  e (R^-S) , 


while  in  place  of  (69) 


- K(R)xV-_-l  - dS  + f k^R)x 


VH(RlR’)dS,  R'  e (R  -S) 


1 

'I 
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The  second  integral  in  each  of  these  two  expressions  can  be 
estimated  immediately.  We  first  note  that 


H (R|  R'  ) 


4tt 


f t , I > 

, ik  R-R' 
1-e  ' 

I I R-R' I 


ik  Y (ik|R-R' | )” 

4it  % (n+1)  1 

n=u 


vh(r1r' ) 


^ (V\R-k'\)  y (n+1)  (ik|R-R’  I)'" 
4tt  ' n=0  ^^+2)  ! 


so  that 


I H ( R I B ' ) I 1 


ke 


k I R-R' 
4it 


l7H(R|R'  ) I < 


4tt 


(87) 


Since  we  are  verifying  the  edge  conditions  the  point  R'  is  a 
point  near  the  boundary  C of  S while  R is  a point  of  S. 
Since  ^ is  bounded,  we  can  draw  a sphere  of  finite  diameter 
D that  contains  S and  R' . Then  |r-R' \ £ D,  and 

, kD  ^ ^ >2  kD 

|H(RrR' ) I 1 ' |VH(R|R'  ) I < . (88) 

With  these  inequalities  we  have  that 


S'- 


, -V  V 


Vq -K(R)  VH  (r1  R' ) +ikZH  (R|  R' ) K (R) 


dS 


Zke 


kD 


— 4 TT 


( lV„*K(R) |+klK(R) I )dS, 


K{R)  xVIi  (R|  R'  )dS 


|K(R)  I I VH(R|R' ) Ids  < 


, 2 kD  (■ 
k e 

4tt  . 


(89) 

]k  (R)  Ids. 

(90) 
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1 


I 

I 

I 

I 


I 

i' 


Since  |k|  and  |Vq*k|  are  integrable  by  assumption,  we  see 
that  the  integrals  involving  H(r|R')  in  (85)  and  (86)  are 
bounded,  so  that  in  verifying  (11)  we  need  concern  ourselves 
only  with  the  remaining  integrals  in  these  two  equations.  That 
integrals  of  this  type  do  satisfy  the  conditions  in  (11)  is  the 
subject  of  Part  I of  this  work. 

To  show  that  the  edge  condition  (iii.b)  is  also  satisfied 
we  use  (69)  and  Theorem  A. 5 to  get  that  at  every  point  r of 
S 

H_(r)-H_^(r)  = -n(r)xK(r).  (91) 

On  any  of  the  curves  C^, 

t(r) • (H_ (r)-H^(?) ) = -K (r ) • t (r) xn (?)  = T(r)-K(r),  (92) 


By  condition  (ii)  in  the  present  theorem  we  see  that  the  limit 
as  e -►  0 exists  and  is  equal  to  zero,  so  that  the  edge  condition 
(iii.b)  is  indeed  satisfied. 

The  boundary  condition  (13)  can  be  shown  to  be  satisfied 
by  first  evaluating  (68)  at  a point  r of  S.  Using  (A. 13)  and 
(A. 14)  we  have  that 


nxE 


= nxE 


S L 


|^VQ-K(R)n(r) xVG (R ] r ) +ikZG (R | r ) n (r ) xK(R)  dS. 


; 1 


I 


( 

t 


\ 

i 

i 


Combining  this  with  (67)  we  obtain  (13) . 


0. 


(96) 


We  now  apply  the  divergence  theorem  to  the  function  E xH  in 


the  region  bounded  by  the  surface  Sq  defined  through  (39)  and 
(40)  (see  also  Figure  1)  and  the  surface  S„  of  a sphere  whose 
radius  R eventually  recedes  to  infinity.  Calling  this  region 
Vq  we  have  that 


(E*xH)dV  = f 


n*  (E  xH)dS. 


V, 


(97) 


I 

I 


The  surface  integrals  over  and  S_^  vanish  in  the  limit 

as  iS  0 because  of  (96)  . The  surface  integral  over 
vanishes  in  the  limit  as  e -*■  0 as  shown  in  Appendix  C.  Thus 
the  integral  over  vanishes  in  the  limit  as  Sq  collapses 

S.  For  the  integral  over  we  use  the  radiation  conditions 

(10)  to  rewrite  the  integrand  as 


to 


n*  (E  xH)  = -E  • (nxH)  = -E  • (Rxh)  = YE*E  + o 


E = YE-E  + O 


IR 

(98) 


since  E = 0(1/R)  (Wilcox,  1956).  Since 


V*(E  XH)  = H*VxE  -E  *VxH  = -ikZH-H  +ikYE>E 


we  have  for  (97) 


ik 


(YE-E*-ZH-H*) dV  = Y 


V 


E-E  dS  + 0 (1) , R 


(99) 


^R 


where  V = lim 
£-<■0 


lim  V 


6->-0 


0 


Since  the  left-hand  side  of  this 


I 


1 

f 

1 


f 

1 


i 

\ 


( 

\ 
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expression  is  purely  imaginary,  we  have  that 

lim  E*E*dS  = 0. 

R->-“  S_ 

K 

For  points  not  on  S,  E satisfies  the  vector  wave 
. 2"^ 

equation  VxVxE-k  E = 0.  From  the  second  of  (8),  definition 
(2),  and  Stokes'  theorem,  we  have  that  V-E  =0  so  that 

V^E+k^E  = 0 . 


Eqs . (100)  and  (101)  imply  that  outside  a sphere  of  radius  R 

enclosing  S,  the  electric  field  vanishes  identically.  This 
is  due  to  the  following  result  by  F.  Rellich  (Muller,  1959). 

Theorem  9:  If  V(R)  is  a solution  of  the  Helmholtz  equation 

V^V+k^V  = 0,  Rek  > 0 
for  R = 1r|  ^ C,  and  if 


fl 

1 V(R) 1 dfi  = o 

U 

w \ 

then  V(R)  vanishes  identically  in  1r|  ^ C 


-►  , 

Thus  E and  H vanish  identically  exterior  to  the  sphere. 
Between  the  sphere  and  S,  they  also  vanish  being  analytic 
solutions  of  (101) . 


(100) 


(101) 
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We  have  then  shown  that  s E2  and  E H2  for  all 

— V 

points  R S.  Since  these  fields  are  continuous  to  the 
surface  S from  and  S_,  these  relations  hold  also 

for  points  R of  S.  We  can  thus  state  the  following  theorem: 

Theorem  8:  The  open  surface  problem  has  at  most  one  solution. 

H . Concluding  remarks. 

The  main  results  of  this  work  are  the  equivalence  and 
uniqueness  theorems  (Theorems  7 and  8,  respectively.)  The 
equivalence  theorem  is  an  immediate  consequence  of  Theorem  5 and 
6,  while  the  proof  of  the  uniqueness  theorem  rests  heavily  upon 
the  results  of  Appendix  C.  In  proving  Theorems  5 and  6 we  note 
that  condition  (iii.b)  in  the  statement  of  the  open  surface  problem 
(Section  C)  was  not  needed  at  all.  What  was  really  needed  was  the 
statement  in  (24)  which  is  a direct  consequence  of  (11).  The 
statement  in  (29),  which  follows  from  (24),  (25),  and  condition 

(iii.b),  is  sufficient  for  the  proofs  but  not  necessary.  Moreover, 
in  deriving  it,  we  only  needed  that  the  function  in  (12. c)  be 
continuous  but  not  Holder-continuous.  The  same  is  true  regarding 
the  existence  of  the  first  partials  of  K in  condition  (v) , 
i.e.  in  proving  Theorems  5 and  6 we  only  needed  the  condition 
that  the  surface  divergence  of  K exists  and  that  its  real  and 
imaginary  parts  are  H-regular  density  functions  for  S. 

The  "excess"  conditions  then  were  necessary  only  in  proving 
the  uniqueness  theorem,  and  specifically  in  showing  that  the 


t : 


I 


{ 
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Appendix  A.  The  behavior  of  certain  integrals. 


In  this  appendix  we  discuss  the  continuity  properties 
of  integrals  of  the  type 


U(R'  ) 


h(R) 
R-R-  1 


dS 


and  V(R') 


g(R)V  — dS, 

J S I R-R'  1 


where  the  function  h is  a regular  density  function  (Part  I, 

Def.  8)  for  S,  while  g is  an  H-regular  density  function 
(Def.  7)  for  S.  It  is  clear  that  U and  V are  infinitely 
differentiable  functions  at  all  points  R'  of  space  which  do 
not  belong  to  S.  The  only  other  points  of  interest  are  those 
points  r which  belong  to  S (but  not  to  its  boundary  C) . The 
proofs  of  the  continuity  properties  of  U and  V at  such  points 
closely  follow  those  for  the  same  type  of  integrals  over  closed 
surfaces.  For  this  reason  we  will  present  here  only  the  results 
and  will  refer  the  reader  to  Muller  (1969)  for  the  proofs. 

Theorem  A . 1 : If  h is  a regular  density  function  for  S,  then 

the  function 


U(R'  ) 


-MiL.  dS, 

S I R-R*  I 


(A.l) 


is  continuous  at  all  points  r e S. 

We  outline  part  of  the  proof  of  this  theorem  to  point  out 

-► 

how  it  differs  from  the  one  for  closed  surfaces.  Let  r e S. 
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j 

{ 


Since  S is  a regular  open  surface,  we  can  erect  a rectangular 
coordinate  system  with  origin  at  r and  the  z-axis  along  the 
normal  n(r)  to  S at  r.  Portion  of  S can  be  described  in 
terms  of  this  coordinate  system  in  the  manner  of  (I.l).  Let 
d be  a positive  number  such  that  the  sphere  |R'-r|  d , and 
with  center  at  r',  contains  only  that  part  of  S that  is 
described  as  in  (I.l)  and,  if  necessary,  we  restrict  d so  that 
no  points  of  the  boundary  C are  contained  in  the  sphere.  We 
denote  this  portion  of  S by  S^(r)  and  let 


U(R') 


h(R) 

•'S-S^(r)  1r-R'  1 


dS  + 


JSd(?) 


dS 

iR-R' 1 


(R'  )+U2  (R'  ) . 

(A. 2) 


We  note  that  exists  as  an  improper  integral  as  far  as  the 

integration  in  the  neighborhood  of  the  boundary  C is  concerned 
because  of  the  properties  of  h.  As  for  the  singularity  at  r, 
Muller  (1969,  lemma  66)  has  shown  that  U^(R')  = 0(d).  Thus 
exists  as  an  improper  integral  if  R'  e S.  For  |R'-r|  £ d/2  and 
lR-r|  £ d,  we  have  that 


1 

[R-R- ( 


R'  -r ' 


I 1 — V — ^ , — 2 

R-R'  R-r  d 


(A. 


so  that 


U- (R’ )-U, (r) I < ^ (R'-rl 
i ■ d 


h(R) 


S-S^(r) 


r 
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This  last  integral  exists  by  virtue  of  the  absolute  integrability 
of  h.  Thus  the  continuity  of  at  r has  been  proven.  The 

proof  for  the  continuity  of  U2  is  given  in  Muller  (1969,  lemma 
69) . Indeed,  the  proofs  in  Muller  and  the  ones  required  here 
differ  mostly  with  respect  to  the  integral  over  S-S^(r).  In 
J 3r  the  integrand  is  continuous,  and  hence  bounded,  while 
here  it  could  become  unbounded  near  the  boundary  C but  it  is 
still  integrable. 

Lemma  A.l:  The  function 

W(R')  = f — dS  (A. 5) 

JS  iR-R'l 


is  continuous  at  every  point  r e S in  the  approach  to  it  from 
either  the  side  or  the  side  S_.  Moreover, 


V — dS 

S+  |M'| 


2irn  (r) 


+ 


dS  + 


dS, 

(A. 6) 


1 

|R-R' I 


dS 


-2-rrn  (r ) 


+ 


dS  + 


dS, 


(A. 7) 


where  by  S^(S_)  we  mean  the  limit  of  the  integral  as  R' 
approaches  r from  the  side  S^(S_),  and  where  Vq  denotes 
the  surface  gradient,  i.e. 


V 


(A. 8) 


1 

c 

'I 

f 
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Proof ; Let  T be  a regular  open  surface  having  the  same  boundary 
C as  S but  no  other  points  in  common  and  such  that  tU  S 
(which  is  a smooth  closed  surface)  has  an  exterior  normal  which 
coincides  with  that  on  S at  their  common  points.  According  to 
Muller  (1969,  Lemma  70),  lemma  A.l  is  true  for  the  closed  surface 
T U S and,  hence,  for  S. 

Remark ; In  lemma  70  of  Muller  it  is  also  proven  that  the  last 
integral  in  (A. 6)  is  continuous  at  all  points  r c (R^-C) . 

Following  Muller  (1969,  Theorems  41  and  42)  we  now  have 
the  following  result 

Theorem  A. 2:  The  function 


V(R')  = f g(R)V  — dS,  (A. 9) 

Js  iR-R’l 


where  g is  an  H-regular  density  function  for  the  surface  S, 
IS  continuous  at  every  point  r e S in  the  approach  to  it  from 


I either  the  side  S,  or  the  side  S_.  Moreover, 
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Moreover,  if  we  let 


H(RlR')  = G(rIr')  + 

4tt|  R-R' 


2 

= - ^[ik  + | + . . . , (A. 11) 


then,  according  to  Muller  (1969,  lemmas  73  and  74),  the  functions 


' r 

h(R)H(R|r)dS  and  g (R) VH (Rj r) dS 

■ s J s 


are  continuous  at  all  points  r of  S,  so  that  together  with 
Theorems  A.l  and  A. 2 we  have  the  following  results. 

Theorem  A. 3 ; If  h is  a regular  density  function  for  the  surface 
S,  then  the  function 


u(R')  = h(R)G(R|R' )dS 
S 


(A. 12) 


is  continuous  at  every  point  r of  S. 


Theorem  A. 4 ; If  g is  an  H-regular  density  function  for  the 
surface  S,  then  the  function 


v(R')  = g (R) VG(R|R' )dS 
•*  S 


(A. 13) 


is  continuous  at  every  point  r of  S in  the  approach  to  it  from 
either  the  side  or  the  side  S_.  Moreover, 


s+(s_) 
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g(R)VG(R|$')dS  = + ^ g(r)n(r)  + 


^ g(?)n  G(R|?)dS  + 
S 


g(r) VpG(R|r)dS  + 


+ I [g(R)-g(r) ]VG(R|r)dS, 
S 


(A. 14) 


where  the  minus  sign  corresponds  to  the  approach  from  S^,  and 
the  plus  from  S_. 

Finally  we  examine  functions  of  the  type 

I g(R) xVG(r|R' )dS, 

J S 

where  the  vector-function  g has  components  which  are  H-regular 
density  functions  for  S.  To  this  end,  and  with  R'  S,  we 


write 


g(R)  XV— ^ — dS  = 
S IR-R'I 


[g  (S) -g  (r)  ] xV— i — dS+g(r)x 

s |r-r'I 


S |r-^' 


dS. 


(A. 15) 


The  first  integral  is  continuous  at  all  points  R'  e (R^-C)  and 
the  proof  is  as  in  Theorem  41  in  Muller  (1969)  . For  the  second 
integral  we  employ  Lemma  A.l  so  that 


g(R)  XV- 


dS  = 


S^(S_)  |R-R' 


(g(R)-g(r)  ]^V— i—  dS±2TTg(r)xn(r)  + 

; lR-^1 


+ g(r) X 


n(R)  ^ 


|R-r| 


dS  + g (r ) X 


Js  R-r 


dS  . 


(A. 16) 
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Using  the  argument  we  employed  for  Theorems  A. 3 and  A. 4,  we  have 

Theorem  A. 5 : If  the  vector-function  g has  components  which 

are  H-regular  density  functions  for  S,  then  the  function 

w(R')  = \ g(R) xVG(R|R' )dS 
JS 

is  continuous  at  every  point  r e S in  the  approach  to  it  from 
either  the  side  S_j_  or  the  side  S_.  Moreover, 


g(R) xVg(r|r' )dS  = ±j  n(r)xg(r)  + 


+ g(r)xv, 


1 

— — dS  + 
l-r\  J 


[g(R) -g(r) ] xV  ■ ^ ^ dS, 
S |R-r| 


where  the  plus  sign  corresponds  to  the  approach  from  S^,  and 
the  minus  from  S . 


r 
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Appendix  B.  An  alternate  derivation  of  the  radiation  conditions. 

In  this  appendix  we  prove  in  a different  way  than  in  Section 
F that  (68)  and  (69)  satisfy  the  radiation  conditions.  For 
R*  ->-00,  we  write 


G(r1r' ) = g(RlR' )+0 


R 


.2 


V'G(rIr')  = ikg(R|R' )R'+0 


,.2 


, (B.l) 


where 


g(R|R' ) 


,ik(R’-R'  *R) 
4tFr’ 


(B.2) 


and  where  (B.l)  are  uniformly  valid  for  all  R e S and  all  R' . 
For  (69)  we  then  have  that 


H(R'  ) = - -^^r-  A + 0 


^ 1 ' 


R' 


, R'  -►  ® 


(B.3) 


with 


A = - 


ik 

4tt 


R' xK(R)e~^^^' ’^dS. 


(B.4) 


For  (68)  we  have,  using  (74), 


E(R')  = - I K(R)  •V'V'G(R|r' )dS+ikZ 


G(r| R' )K(R)dS.  (B.5) 


But,  for  R'  -+■ 


1 

1 5 ^ ^ 

■ 1 ■ 

'7'V'G(r1R')  = ikV’g(R|  R' ) R'+O 

= -k  g(R|R' )R'R'+0 

R'  ^ 

so  that 


E(R') 


ikZ 


R’ (R' •K)gdS+ikZ 
S 


( gKdS+0 

f \ 

1 

, 2 

Js 

^ J 

= -Z 


[ [K(R) xR' ]gdS+0 

1 

“77 

1 

[r- 

, 1 

A. 

1 

— R'xA+0 

R'2 

k J 

9 

R'  -*■ 

From  (B.3)  and  (B.7)  the  radiation  conditions  (10)  follow 
immediately. 
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Appendix  C.  The  estimation  of  an  integral. 

In  this  appendix  we  will  estimate  the  intjsgral 

n*(S*xH)dS,  (C.l) 

with  defined  in  (40) , which  appears-  in  Section  G.  To  this 

end  we  choose  a point  Mq  of  C and  with  it  as  origin  we  erect 
a rectangular  coordinate  system  with  the  x-axis  in  the  direction 
of  the  tangent  vector  to  C at  Mq,  and  the  z-axis  in  the 
direction  of  the  normal  to  S at  M.  The  y-axis  then  lies  on 
‘the  plane  tangent  to  ^ at  M and  is  directed  toward  S.  On 
the  yz-plane  we  erect  a polar  coordinate  system  (p '#<!>')  as 
described  in  Section  B.  The  angle  (|>'  is  measured  coxanter- 
clockwise  from  the  positive  y-axis.  On  this  plane  we  draw  the-; 
circle  p'  e > 0 and  we  estimate  the  integrand  of,. (C.l)  at 
points  of  this  circle  except  at  the  point  which  belongs  to  the 
surface. 

it  . '' 

Since  the  unit  normal  n is  pointed  toward  M,  we  have  that 
n*  (E*xH)  = -p'*(E  xH)  = -(y  cos  <f»'+z  sin  4)')*(E  xH)  = 

= - cos  (()' (£*x5)y-sin  4>'(E*xh)^  = 

= -cos  <T(E/H^-E/H^)-sin  * (E^Hy-E^X^  = 

= E„*  (cos  (J)'H  -sin  (j)'H  )+H  (sin  <f’E  *-cos  (^'E  *) 

X 2 Jr  ^ J “ 


(C.2) 
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As  mentioned  above,  the  integral  in  (C.l)  appears  in  Section  G 
where  the  functions  E and  H are  assumed  to  be  solutions  of 
the  open  surface  problem.  Then,  by  the  equivalence  theorem  7, 

E and  H are  given  by  (68)  and  (69)  , respectively.  These 
expressions  appear  rewritten  in  (85)  and  (86)  in  a way  so  as 
to  isolate  the  singularity  of  G.  As  pointed  out  in  the  disucssion 
there,  in  examining  the  behavior  of  E and  H near  the  edge 
we  need  only  cosider  the  first  integrals  in  (85)  and  (86).  Thus, 
with  R'  = (p ',(()'),  we  have 


E(R')  = 


- U: 


if 


ikZ 


R-R' 


R-R' 


K(R) 


dS  + 0(1) 


(C.3) 


H(R')  = - 


4tt 


K(R)  XV 


R-R' 


dS  + 0 (1) 


(C.4) 


We  will  estimate  (C.2)  by  using  these  last  two  expressions.  We 
note  that  the  densities  appearing  in  them  are  H-regular  density 
functions  for  S by  condition  (i)  in  Theorem  6 so  that  we  can 
use  the  results  in  Part  I.  Though  K and  are  complex-valued 

we  will  treat  them  as  real-valued,  the  extension  to  complex- 
values  being  obvious.  From  Theorem  1,  Part  I,  we  have  that 


E^(R')  = 0(1) , 


p’  0. 


(C.5) 


We  will  now  prove  that  the  same  is  true  for  H^.  From  (C.4)  and 
with  respect  to  the  coordinate  system  above,  we  write 
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- k Jg  h 


- K (R)  dS  + 0(1) 


[K  (R-)(z-z')-K  (S)(y-Y*)]  -f,  -3 
^ J S ^ ^ R-R' 


+ 0(1) 


(C.6) 


We  note  that  the  integral  here  is  a combination  of  the  integrals 
appearing  in  (1.43)  and  (1.74)  and  will  be  treated  in  the  same 
way  up  to  a certain  point.  As  remarked  at  the  beginning  of 
Section  D,  Part  I,  the  integration  over  S can  be  split  into 
two  parts  one  over  S(A)  and  one  over  S-S(A),  where  S(A)  is 
that  neighborhood  of  S about  Mq  whose  projection  on  the 
xy-plane  is  the  region  A defined  in  (1.12).  The  integral 
over  S-S(A)  is  continuous  in  a neighborhood  of  Mq,  and  its 
limit  as  p'  ^ 0 is  equal  to  the  analogous  integral  obtained 
when  R is  replaced  by  5.  The  integral  over  S(A)  can  be 


written  as 


[K  (R) (y-y')-K  (R) (z-z') ] — ^ 
S(A)  ^ y I R-R 


dS 

I 3 


Kj,  (y-y ' ) sec  <J^(x,y) 

I R-R' I ^ 


dxdy  - 


K^(x,y)  (y-y')sec  ip(x,y) 


I R-R' 


dxdy  = 


(C.7) 


where  ij;  is  defined  in  (1.25)  and  is  the  angle  the  z-axis  makes 


with  the  normal  to  the  surface  at  (x,y) . 


I 
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For  we  write 


[K  (x,y)sec  iMx,y)-K  (0,y)sec  ij^(0,y)]  dxdy  + 


R-R'  ' ^ 


k : 

r i 

(•  i 

f 

I 

' I 

• i 


\ 


' 

K„(0,y)sec  i|'{0,y)  (y-y') 

1 1 
1 , 3 3 

Ja 

J R-R'  1 R _ 

dxdy  + 


1 K (0,y)sec  <jj(0,y)  dxdy, 

J A ^ R-^ 


(C.8) 


where 


R =Vx^+(y-y’)^+(F(0,y)-z')^  . (C.9) 

The  second  integral  above  has  been  discussed  in  Part  I,  eqs. 

( 1 . 32 ) - (1 . 36 ) , and  has  been  found  to  be  bounded.  For  the  first 
integral,  call  it  I^,  we  write 


^3  = 


[K  (x,y) -K  (0  ,y)  ] sec  'J^(x,y)  - dxdy  + 

A ^ ^ Ir-r'I 


[sec  \p(x,y)-sec  (0,y)  ]K  (0,y)  ^ ^ dxdy  . 

A I R-R' I 

The  second  integral  here  can  be  shown  to  be  bounded  by  treating 


(C.IO) 


(C.ll) 
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Since  9K^/3x  is  an  H-regular  density  function  for  S,  the 
integral  reduces  to  one  as  in  (1.29)  which  has  been  shown  to 
be  bounded. 

For  I2  we  write 


^2  = 


Z-2  ' 


tK  (x,y)sec  ip(x,y)-K  (0,y)sec  i|j(0,y)]  — - dxdy  + 

y y Ir_r'  I 


K„(0,y)sec  (0  ,y)  (z-z  ' ) 


A 


R-R'l^  R^ 


dxdy  + 


K,,(0,y)sec  ijj(0,y)  ^ ^ >7^  dxdy  + 


K (0,y)sec  i|;(0,y)  dxdy. 


(C 


R' 


The  first  two  integrals  in  this  expression  can  be  treated  as  the 
corresponding  ones  for  . The  third  integral  is  of  the  same 
type  as  the  one  appearing  in  (1.76)  and  is  bounded.  We  can  then 
write 


^1  ^2 


[K^(0,y)  (y-y'  )-K  (0,y)  (f  (O  ,y) -z  ’ ) ] dxdy  + 0(1), 


R' 


p'-*-  0, 


(C 


.12) 


.13) 


In  order  to  estimate  this  integral  we  will  first  rewrite  it  slightly 
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differently.  We  note  that  K and  K are  evaluated  at 

z y 

points  of  S which  also  belong  to  the  yz-plane,  i.e.  on  the 
curve  z = F(0,y).  As  explained  in  Part  I,  after  eg.  (1.2), 
the  normal  to  the  surface  at  such  points  and  the  normal  to  the 
curve  coincide  because  of  condition  (1.2).  This  normal  is  given  by 


-F  (Ofy)y+z 


a+Fy^(0,y) 


(C.14) 


The  unit  tangent  vector  to  the  curve  is 


y+F  (0 ,y) z 


(C.15) 


a+F^"-(0,y) 


so  that  Txn  = x.  Since  the  current  density  K is  a surface  field, 
it  can  be  written  as 


K(0,y)  = K,j,(0,y)T+K^(0,y)x, 


(C.16) 


so  that 


KT,(0,y) 

Ky(0,y)  = - ■: 

/l+F^^(0,y) 


K^(0,y)  = 


F (0,y)K^(0,y) 


/l+Fy^(0,y) 


(C.17) 


and,  since  sec  i|^(0,y)  = /l+F^  (0,y),  we  have  for  (C.13) 


K (0,y) 

= [F  (0,y)  (y-y')-(F(0,y)-z')  ] -!^-3 dxdy+0(l),  p ’ 0 


We  now  employ  the  coordinates  (s,p)  and  use  the  results  for 
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the  integrals  in  (1.45)  and  (1.78)  to  get 


I 


1 


•b  K^(0,p)  (1-Kp) 

0 (p-y')^+(F(0,p)-z')^ 


[Fp (0, p) (p-y ' )-(F(0,p)-z') ]dp+0 (1) , 

(C.19) 


where  k = k(0),  k being  the  curvature  of  the  projection  of 
the  boundary  C on  the  xy-plane  defined  in  (1. 11).  With 


u = p + i(F(0,p),  w = y'+iz'  = p'(cos  (})'+i  sin  0')  (C.20) 


we  can  write 


I^-l2  = 2Im 


K^(0,p) (1-Kp] 
u-w 


du  + 0(1) 


(C.21) 


where  is  the  curve  whose  equation  is  z = F(0,p), 

0 p _<  b.  Since  K^p  is  Holder-continuous  on  C^,  including 
the  end-points,  and  since  K^(0,0)  = 0,  we  have  that  the  integral 
is  a bounded  function  of  w tending  to  a definite  limit  as  w 
approaches  the  end-point  (0,0)  along  any  path  (Mus)chelishvili , 
1953).  From  (C.21),  (C.7),  and  (C.6)  we  have  that 


H^(R')  = 0(1) , p'  ^ 0, 
so  that  (C.2)  becomes 


(C.22) 


|n-(E*xH)|  < a(|h^|  + |h^|  + |e^1  + |e  1)  < a(|h|  + 1e|),  p'  -V  0,  (C.23) 

y z y ^ 


where  A is  a positive  real  number. 


From  the  Cauchy-Schwarz  inequality  and  (C.23),  we  have  for 


(C.l) 


I "■ 


(E*xH)dS 


2 r 
< 


I 1 

n*  (E  xH  dS 


2tt 


->-*  I 

n*E  xH  e 


1 I dt 


d(}) ' ds ' 


< e 


' rL 

f2Tr 

. O' 

0 ■ 

1-ep 


. 

ds ' 


d(f> ' ds ' 


2it 


^ ^ 2 1 
n *E  XH  I d(}) ' ds  ' > < 


2 

< A'e 


27T 


( lEl2+lHl^)d(t)’ds'  , 


where  A'  is  a positive  real  number.  Combining  this  result 
with  (11)  , we  get 

n*  (E*xH)dS  = 0 (e")  , e 0. 

E 


(C.24) 


(C.25) 
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